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Abstract 

The two-loop massive operator matrix elements for the fermionic local twist-2 operators 
with external massive fermion lines in Quantum Electrodynamics (QED) are calculated 
up to the constant terms in the dimensional parameter e = D — A. We investigate the 
hypothesis of Ref. pQ that the 2-loop QED initial state corrections to e^e~ annihilation 
into a virtual neutral gauge boson, except power corrections of 0((m^/s) ), A; > 1, can be 
represented in terms of these matrix elements and the massless 2-loop Wilson coefficients 
of the Drell-Yan process. 



Alexander- von- Humboldt Fellow, 
^deceased. 



1 Introduction 

The QED corrections for differential distributions in e^e~ annihilation and other high energy 
reactions in which electrons or positrons participate, are particularly large due to the presence of 
physical logarithms ln(M^/r7ig), with M a characteristic scale of the process and nie the electron 
massjfl Therefore, it is necessary to account for the QED initial state corrections up to O(a^) for 
precision measurements in the various energy regimes in e^e~ annihilation having been explored 
so far and those which are planned to be investigated in the future, cf. [3HT2]. For the corrections 
to the inclusive Born cross section cr(s), with s the center of mass (cms) energy squared, power 
corrections oc {m1/s)^,k > 1 can be safely disregarded. While the 0{a) corrections are known 
for a large number of reactions, cf. [13], the corrections beyond the universal contributions 
0((aL)'^), 1 < fc < 5 [HllIS], to higher orders, were only calculated analytically once at 2-loop 
order in Ref. pTj.u Besides the logarithmic orders 0{a'^L'^, a^L) with L = ln(s/mg), the constant 
terms Ola"^) are of interest. 

In Ref. [Ij it has been proposed that the 2-loop corrections can be calculated using a 
factorization-representation of the scattering cross section in terms of Mellin convolutions of 
massive local operator matrix elements (OMEs) and the corresponding massless Wilson coeffi- 
cients, which are those of the massless Drell-Yan process [20l[2T] . The operator matrix elements 
are formed by the local twist-2 fermionic operators, being obtained in the light cone formal- 
ism [22] , between massive on-shell electron states. These operator matrix elements bear all the 
mass dependence, fi'^/ml, and are universal quantities. The massless Wilson coefficients account 
for the process dependence and are functions of the ratio s//i^, with s the sub-system cms energy 
squared. Here, /x^ denotes a factorization scale. 

In case of the heavy flavor corrections to deep-inelastic scattering the above method was used 
to calculate the massive Wilson coefficients in the region Q^ ^ 771%, with Q^ the virtuality of the 
exchanged gauge boson and rng the heavy quark mass. It has been shown that this description 
yields all but the power suppressed contributions in Refs. [THl[19] comparing to the complete 
semi-analytic calculation [23] at O(a^) U- In this case the massive OMEs are formed between 
massless on-shell quark and gluon states. In Refs. [261428] 2- and 3-loop heavy flavor corrections 
for different unpolarized and polarized nucleon structure functions and transversity, respectively 
for their moments, have been calculated. 

In the present paper we compute the O(a^) local OMEs with massive external fermions 
in QED up to the constant part. As a by-product we obtain the QED contributions to the 
2-loop non-singlet and pure-singlet anomalous dimensions, within a massive calculation. We 
investigate, to which extent the decomposition, having been proposed in Ref. [1], in terms of 
massive local OMEs and massless Wilson coefficients is possible. 

The paper is organized as follows. In Section 2 we summarize the decomposition of the initial 
state corrections to the inclusive e~^e~ annihilation cross section into neutral vector bosons to 
0{a'^) using the renormalization group method, cf. [Ij. The renormalization of the OMEs is 
described in Section 3. In Section 4 we compute the OMEs to 0{a). The details of the calculation 
of the O(a^) corrections to the OMEs are given in Section 5. We discuss the structure of the 
contributions to the differential scattering cross section dae+e-/ds and compare to the results 
in the literature. Section 6 contains the conclusions. Technical aspects are summarized in the 



^ This also applies to the QED corrections in eA^ scattering, cf. [2]. 

There are only a few other complete analytic 0{a^) two-loop calculations with massive particles available. 
Examples are the heavy quark fragmentation functions [TB], the /i-lepton decay spectrum [T7] to 0{a^L), while the 
O(a^) contribution was given numerically, and the heavy flavor Wilson coefficients for large virtualities [1^119) . 

^For a fast and precise numerical implementations of these corrections in Mellin space cf. '24' . 



appendices. 



2 The Renormalization Group Method 

The O(a^) QED initial state corrections to the e^e~ annihilation cross section into virtual neutral 
gauge bosons (7*, Z*) in the limit s 3> ml can be expressed in the following form [1] : 
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ds' ds' ds' ds' 

Here, s' denotes the invariant mass of the virtual vector boson and s the cms energy of the 
process, 



s' = xs, X G [0, 1] . 



(2) 



The term I refers to the photon radiation contributions, II to the flavor non-singlet contribution 
due to fermion pair production, and III to the corresponding flavor pure-singlet contribution. |^ 
One may represent the scattering cross section in Mellin space by applying the integral transform 
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Since in the present calculation power corrections of 0{ml/s) are disregarded, the following 
principle structure [1] with respect to the scales ml, s and a factorization scale yU^ is obtained, [1] : 
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To 0(a;^) the last process in ([6]) does not contribute. Here ao{s) denotes the Born cross section 
for e^e~ annihilation into a virtual gauge boson (7, Z) which decays into a fermion pair //, see 
e.g. |29], 
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^In Ref. [T] four contributions were considered dividing those to process I into two pieces according to the 
genuine 2 — >■ 3 particle scattering cross sections. 
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Here a denotes the fine structure constant, with a = Ana, Nqj is the number of colors of the 
final state fermion, with Ncj = 1 for colorless fermions, s is the cms energy, fi is the spherical 
angle, 9 the cms scattering angle, and the effective couplings G'j(s)|j=i...3 read 
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The reduced Z-propagator is given by 
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where Mz and Tz are the mass and the with of the Z-boson and mj is the mass of the final 
state fermion. Qej are the electromagnetic charges of the electron {Qe = —1) and the final state 
fermion, resp., and the electroweak couplings Vi and a^ read 
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where 9^ is the weak mixing angle, and /^ ^ = ±1/2 the third component of the weak isospin for 
up and down particles, respectively. 

The factorization mass cancels in the physical cross section in each order of the coupling 
constant. The initial state fermion mass dependence is solely encoded in Tu. The operator 
matrix elements are given by 
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where Op ' and Oy are the local twist-2 fermion and photon operators, 
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for the fermionic non-singlet (NS), singlet (S), and photonic case, [30]. Here, S is the sym- 
metrization operator of the Lorentz indices fii, . . . , fiN, ip denotes the electron field, F^j/ the 



photon field-strength tensor, and D^ = d^ — ieA^ the covariant derivative, with e = ^^/{4:7^)'^a 
the electric charge and A^ the 4-potential. We consider only one fermion species. For ( !T7|) both 
the flavor non-singlet (NS) and pure singlet (PS) terms contribute. It turns out that also in 
the present case the contributing functions both for the massive OMEs and the massless Wilson 
coefficients |31j can be related to nested harmonic sums [32ll33] and multiple zeta values [Mj . 
These structures simplify further applying algebraic and structural relations [35] . One may per- 
form the calculation in Mellin space completely and represent the result in x-space analytically 
resp. numerically performing the inverse Mellin transformation [35|l36]. 
The following representations apply 
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Here i,j denote the external particles in the scattering process. The coefficients a„r- and bnr of 
the above series adjust such, that the physical cross section is independent of /i. 

r^j, aik and the scattering cross sections (Tjj obey the renormalization group equations 
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The /3-function f3{g) is given by 
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with g the electromagnetic coupling. In QED the first expansion coefficients are [38] 
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in the case of one light fermion. The running coupling alfi"^) 
is obtained as the solution of 
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with ao = a{m.l). We rewrite the renormalization group equations replacing jid/dn by 2d/dL, 
resp. —2d/dX, with 
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where a = a{fi^). Here the anomalous dimensions 7ij(A^, a) have the series expansion 
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For later convenience we also introduce the splitting functions in A^-space, 
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The higher expansion coefficients of the OMEs, T'^j, T^j and Tjj, ( !79| I142H144I) . are defined such 
that they do not flip sign under the Mellin transform. 

For the process under consideration one obtains to Ola"^) : 



Fee ( N,a,^ 
mi 



l + a 



+a^ 



-^^iN)L^V^iN) 



\lfXN) {^fXN) - 2/3o) + ^SH^)7^?(iV) j L 



+\ {-7£HiV) + 2/3or(°) - lfKN)T^^){N) - ^f:iTfAN)]L 






+ Oia") 



(36) 



ae.e. [ N,a,— 



l + a[-^?XN)X + afXN)\ 



+a^ 



\lfKN) (7(?(iV) + /3o) + \lf:j{N)^^^KN) ) \' 



+ {-lil\N) - /3oa(°) - 7l?(iV)a(°)(iV) - ^\l^a^^{N)} X 

+ Oia'') 
+ 0(a^) , 







+ 0(a3), 


F^e AT, a, ^ 


= a 


-Ii^'Kn)l + t(^} 


/ s'\ 


= a 


1 

_~2 


7SHiV)A + a(;) 



(37) 

(38) 
(39) 



with the 1- and 2-loop splitting functions given in Refs. [391 SO]. We express the coupling 
constant a = a{^^) by (!30!) and assemble the differential scattering cross section ([1]) in terms of 
the three contributions: the flavor non-singlet terms with a single fermion line (I), those with an 



additional closed fermion line (II), and the pure-singlet terms (III). Here we corrected Eq. (4.19, 
4.22, 4.26) in Ref. \l\ 0. 
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We will later rewrite Eqs. fH0l - |42|) in terms of L. 
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3 The Renormalization of the Operator Matrix Elements 

Since the massless sub-system scattering cross sections aab{,z,s/ n'^) to O^a^) are known, the 
corresponding massive operator matrix elements for e^ — )■ e^,e^ ~^ 7 and 7 — )■ e^ transitions 
have to be calculated. The latter two processes contribute only to first order of the coupling 
constant. Here, the external fermion is a massive particle, contrary to the cases studied in [261428] . 
The bare OMEs are given by 
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with a the unrenormalized coupling constant. The computation is performed in D = 4 + 6 
dimensions. The electron mass is renormalized on-shell 



p^ = ml{fi = rrie) 
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with p the momentum of the external fermion. Thus renormalization concerns the wave function, 
the charge renormalization, and the ultraviolet singularities of the local operators. For process I 
counter terms emerge at 0{d?). Due to the finite fermion mass, no collinear singularities emerge. 

i) Wave function renormalization. 

The bare wave function ipQ is renormalized by 

i^o = v/^MV'- (46) 



^The contributions to process II have been also calculated in Ref. [25j . 



The Z-factors are obtained from the fermion self-energy, see Figure 1, 
Expanding ^ [p, rrie) around p"^ = ml one obtains 



Y, = ^ + ^rnl—E,{p\m^ 



+ ^2(p',me)|^,^^. 



(47) 



(48) 






(a) 



(b) (c) 

Figure 1: The self-energy diagrams. 
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To O(a^) Z2 is given by ^41] 
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with the spherical factor 
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and (k denotes Riemann's (^-function at integer values k > 2. In the MS scheme the factors 
Ss are set to one at the end of the calculation. In ( H9l) we separated the terms contributing to 
processes I and II. At 2-loop order also counter terms (CT) contribute with the Z-factor Zqt to 
the OME Agg, which will be calculated in Section [51 

After wave function renormalization and accounting for counter terms the OME is denoted 
by Aee- Up to O(d^) it is given by : 
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ii) Charge renormalization. 

The bare coupling a and the renormahzed couphng in the MS scheme are related by 
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/3o(^/) = -3^/ • (54) 

The above relations would apply to nj manifestly massless fermions. Since the fermion lines are 
all massive the coupling constant is first being obtained in a MOM-scheme, which is defined by 
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cf. [27]. Here II h,bf {0,1^1) denotes the on-shell vacuum polarisation calculated using the 
background- field method [12], 
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Finally, we transform back to the MS scheme using 
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iii) Renormalization of the composite operators. 

We express the inverse Z-factors in the MOM-scheme : 
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The transformation to the MS-scheme is obtained by fl62|) . We spht the OME into the parts 
I-III. The corresponding Z-factors are given by 
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The OMEs to two-loop order after wave function and charge renormahzation are given by 
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The renormahzed OMEs in the MS-scheme are finally given by 
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Here T-- and F^ fl90] - [93|) denote the constant and linear term in e of the unrenormahzed one- 
loop massive OMEs, and T\j (I142H144|) the corresponding constant part of the two- loop OMEs. 
Also the OMEs A^^yiN) and A^e{N) have to be calculated to 0{a). Here only the operators 
have to be renormalized. 
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4 The 0(a) Operator Matrix Elements 

The massive operator matrix element for process I, Aie , emerges at 0(a) and is given by f rM|) . 
In the calculation of process HI to 0{a^) also the OMEs Ai-y and A\J , (175| [76]) contribute. The 
Feynman rules for the operator insertions are given in Figure [21 cf. [27]. The external lines are 
taken on-shell, i.e. p"^ = ml for the fermion and p^ = for the photon lines, and the vector A is 
light-like with A ■ A = 0. 

In the following we present the OMEs in x-space. The unrenormahzed OMEs are given by 
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where for AiJ (x^e) the contributions up to 0{e) are needed. In Figures [3H1] the diagrams 
are shown, except self-energy diagrams, which contribute to the 0{a) OMEs, cf. Figure 1. 
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The expansion coefficients are the leading order splitting functions P^^- [39j, F^y, and T^j , 
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Figure 2: Feynman rules. 



We define 
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The +-prescription, used to regularize some of the terms, reads 
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and g{x) e T>[0, 1] denotes a test function, [43j. 

The sphtting functions obey the well-known relations 

Plfix) = Pi'iKl-x), X<1 



-'67 i"^/' ~ -' 67 \^ 



X) 



and 



rfxPW'^s-fx) 



rfxx[Pi°)(x)+pW(x)] 



0, V/ G N 
. 



Eq. (1871) derives from fermion number conservation with 

Eq. flHHl) results from the conservation of 4-momentum. 
The 0{e^) terms are 

rW(x) = -8T>i(a;)-4Do(a;)+4(5(l-a;) + 2(l + x)[21n(l-x) + l] 
'1+x^ ( , , 1 

-4 Y-;\Hi-^) + -, 










-2 ^^^^ ^^' [21n(x) + l] 



(83) 
(84) 

(85) 
(86) 

(87) 



^9) 



X 



(90) 
(91) 
(92) 
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Figure 3: Diagrams contributing to A, 



(1) 

ee ■ 





Figure 4: Diagrams contributing to A\J and A, 



+ crossed 



ie7 ■ 



.(0) 



T^(0), 



for re7 (x) cf. [T8l[T9]. Here we corrected Eq. (4.28c) in \T\. The linear term in e, T^^ (x), reads 



rl°V) 



-{ 42)2(x) + 42)i(a;) + C22?o(a;) + ( 4 + -C2 ) 5(1 - x) 



-2(1 + a;; 
'l + x^ 



-2 



X 



ln^(l-a;) + ln(l - x) + -(2 



ln2(l - x) + ln(l - x) + -C2 



(93) 



In the differential cross sections fH0] - [42l) different convolutions of the expansion coefficients of 
the leading order OMEs and the leading order Drell-Yan scattering cross sections occur, see also 
Appendix A. Unlike the expansion coefficients of the massive OMEs, the coefficient functions 



(jj- are process-dependent quantities. In case of the e"'"e annihilation process the massless 



u 



Wilson coefficients can be obtained from those of the QCD Drell-Yan process given in [20] 
to 0{a?)^ adjusting the color factors. The 0{a) Wilson coefficients read a 



^1?(^) 



^e^a:) 



16Di(x) - 8(2 - (2)^(1 -x)- 8-^^ - 4(1 + x) [2 ln(l - x) - ln(x)] 

-L Jb 

-P^){x) [2 ln(l - x) - ln(x)] + 1 + 6x - Tx^ . 



(94) 



The combination 



2ri°Hx) + a(°)(x) 



-82)o(x) + 8(C2-l)<5(l-a;)-8 



ln(x) 
1 — X 



+ 4(l + x)(ln(x) + 1) 



_4i±^ [1 + ln(x)] - 8(1 - C2)'5(l - x) 
1 — X 



(95) 



occurs in (140] - [42|) . Here the logarithm L still bears a x-dependence. Referring to L instead, the 
0(ao) contribution to the annihilation cross section (l40l) is given by 



■^1,1 






(96) 



^Note a typographical error in Table 6, Ref. [44] 
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which resembles the well-known behaviour of the splitting function Pee (x) contributing through 



5 The 0{a^) operator matrix elements 

In the following we discuss the Ola^) contributions to operator matrix elements for the processes 
I-lII individually. Further, we investigate their contribution to the differential scattering cross 
sections (|10ti42D. 



5.1 The OME A^^J^^ 



ee 



rn\ T 

The Feynman diagrams contributing to Aee are shown in Figure except those contributing 
to the wave function renormalization, cf. Section 3. Furthermore, the counter terms shown in 




(f) (S) (h) (i) (J) 

Figure 5: Feynman diagrams for the calculation of the massive two-loop operator matrix elements 

.(2),I 

Figure |6] contribute. The corresponding Z-factor is given by 

.2\e , j2_, ^ 24 



^ct(^) 



Tfi 

-i ] \ -^-^(l - ^) - — 25o(a;) - 24Di(a;) + 162)o(x) 

/i / I £ 6 



- (64 + I8C2) 5(1 - s) - 12 - iV 



24 
8 ) (5(1 - x) - 24Do(x) 



■ (97) 




Figure 6: Counterterm diagrams. The black squares represent the counterterm vertices, cf. j45j . 



^ (21 I 

It contains a term oc A^, which cancels a corresponding contribution in Aee ■ It can be 



rewritten, cf. 



usmg 

dxN x^^^5{l - x) 



j dx (-^^A 5{l-x) = ~ j dxx^~' [x6'{l - x)] 



(9^ 
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and similar relations. 

The two-loop diagrams shown in Figure [5] are calculated using the Feynman rules given in 
Figure |2l For example, the ladder diagram in Figure |5]d yields : 

J (27r)^(27r)^ D'^W^W^W^W^' I " 2j , {^^) 

where n = N — 1. We introduce the following short-hand notation for the denominators : 

Di = kl-m\ D2 = kl-m\ D3 = (fci - p^ , 

D4 = (ki - k^f , D5 = {h -ki+ pf - m" . (100) 

The same propagator structure can be obtained for all diagrams from E^ to Eh by choosing the 
flow of momenta appropriately, while conveniently keeping all of the dot products coming from 
the Feynman rules in Figure [2] as simple as possible, i.e., as powers of A • fci, A ■ ^2 and A -p. The 
remaining diagrams are calculated directly. The Dirac-structure of the numerators is projected 
multiplying by 

\{^ + m) (101) 



and taking the trace. We applied FORM |16] for this calculation. This produces a linear combi- 
nation of products of all possible dot products of A, fci, k2 and p. After canceling as much as 
possible these dot products against the propagators, and choosing ^2 -p as the only remaining ir- 
reducible numerator not involving A, diagrams [5K to[5li can be expressed as linear combinations 
of the following type of integrals 

d^ki d^k2 (A ■ A;i)'^(A ■ k2 
^uiU2U3U4,us ^ I (n„\D (o^\D n'^i rt^2 pi'^B n'^i n'^5 ' 



Aa,b ^ ; ^ ■■^1 ^ ■■^A v— •"!/ V" -z/ ('1021 



a,b ^ / "- "^1 »- A.2 n.2-pK^-i^l) K^ ' i^2) /-,^o\ 

^1^1^2^31/4^5 I (n \D (OTr\D T^'^1 T^^-i T^i^B T^i^i n^^r, ■ l-'-'-''^J 



d'^ki d^k2 k2 ■ p{A ■ ki)%A ■ k2 

Diagrams [5^ and O can be written as linear combinations of integrals with the structure 

F-^^ = f ^""^^ ^""^^ (^ • ^^^"(^ • ^'^' VfA kV(A kr-'-^ no4i 

U1U2U3U4U, J {2tx)D{2ti)d Dl^D^^^DfDl'Dl^ ^^ iM 2; , K ) 

and diagrams ^g and |5]i are given by linear combinations of integrals of the form 



^b "-1 



F-' = I ^""^^ ^""^^ (A • A:0"(A ■ fc2) ^ , 

^ VxV2V-iV4VT, / (O^^D (r)^\D n'^in'^2 pjl/ip,Uip,U5 /_^\^ fJ V^ '^U ■ K-^^'^J 



3=0 

We list now the results for diagrams |5^ to [5)i in terms of the integrals above : 
Diagram [S^: 

2 (-^ ~ 2)(-^ ~ ^) [^oiiii + ^liioi ~ ^liiio ~ ^lioii + 2^02110 ~ ^A 



0,1+n 
12001 



iO,l+n] 
M2101 I 



- (A • P) A°i1,,] - 2{D - A)ml [(A ■ p) A?f,,, - A?i\tr] + 8m^ [A^^^^ + A\ 

+ \{D -2){D- 8) [A?il,tr - <\ti" + (A ■ ;.) «ni - ^"loi)] 

~2(-D — 2) [Aj^^j^oi ~ (A • p) ^lioii ~ ^liiio + 2-Bi2ioi + 2-Bi2oii J 

-^rnl [A°^\|r + ^?2\to ] - 16m^A?,\tr (106) 
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Diagram [Sb: 



Diagram Efc: 



Diagram Eli: 



l^i^ — ZJ [/liiiio — ^11110 ~ {^ ' P) ^11110 ~ ^02110 + ^-Dl2110 J 

+4(D - 2)ml [A°i\to - ^2ino] - ^D' - AD + 8)mlA%\\-, - 16mtA%\r, (107) 
{D -2){D- 8) [(A • p) A-f,,, - (A ■ p) <,o, + A^^^^ + Aj+lf - ^omf 

^11011 + ^OllllJ + ^\^ ^J^llllO + °^e [^^21011 + ^^21101 ^21110 

-^2oYif] - (/^ - 2)(D - 4) [2A^+;f + {A.p) <,,,] - 32m^A^+';f 
+4(D — 2) [A;^oiii + ^liioi + (^ ■ P) (^lioii ~ ^liiio) ~ ^20101 ~ ^liiio 
-2Bltrof - 2Blt;n'] + A{D - A)ml [2A}+';f - {A ■ p) <,,, - A^^il.J (108) 



[D — 2) [2yl2iio'o — "^A^-^^^^Q + A]^{]^;^Q + (A ■ p) (Aj^^j^^Q — Agiioo) + 2-B21110 J 
-2{D - 4)2m,2A^riTo - ^{D - 2)ml [(A ■ p) <no + 2^^+^ + <iio] 
-32m^A^+To (109) 



Diagram Efe: 



2{D - 2) (A • p)^ [E,7,i, - eZ,, + EoTiii] - 2D (A • p) E^,, + 4E?foii 
+2(D - 2) (A • p) [Elf,,, + EiYoii - E\f,,] + 2(D - 6) [eH,, - EoYiii] 
+2(D - 4) Kini - i^iioii + ^mio " ^moi + (A " P) <ni] " 4i5,°'' 



Diagram OF: 



+4(D - 3) (A • p) [E\f,, - i^o^lii] + 4 (A ■ p) [Eli,, ' ^?iioi] + ^E^,, 
-4El',\„ + 8ml [{A . p) {E',[\„ + Elf„,) + eI[\„ - E^f,,,] (110) 

AiD - 2) [{A . p) E',i\„ - Eli\„] + l6mlEl{\„ (111) 



Diagram [S^: 



2(D - 4) [(A ■ p) F^l,, - F^il, - F^i%, + F^f„, + F^',,,] + 2D (A ■ p) F^° 



10111 



+2[U 2) [[A ■ p) [r ,,,Q, -T 10111 j + -^11101 -f'liiio + -^oiiii -^lioiij 

+A{D - 3) [F^il,, - F^l,,] + 2{D - 6) (A • p) [F^f,,, - F^',,,] 

+4 (A . p)^ [F^i%, - FZn - i^o°i?ii] + 4 (A ■ ;.) [F^f„, - F^f,,, - F',il„] 

+8 (A ■ p) ml [Flf„, + F',i\„] - 8mlF^f„, + 8mlF^il„ (112) 

Diagram [5]i: 

4(D - 2) (A • p) [F^i\„ - F^f„, + F^i%,] + 16 (A ■ p) mlF^f„, . (113) 

Various of the integrals appearing in these expressions have only three or four propagators. 
The 4-propagator integrals can be represented in terms of up to three Feynman parameter 
integrals over the unit cube. In some cases, a direct calculation will give integrals with the 
following structure 

I{e,n) = / dx dy dz x'^ f {x , y , z; e) , (114) 



^0 ^0 
15 



while in other cases they will be of the form 

I{e,n) = / dx dy dz x'^y'^f{x,y, z;e) . (115) 

Jo Jo Jo 

In the first case, the integrals represent a Mellin transform, and only the y and z integrals have 
to be performed. A mapping of Feynman parameters providing this case can be applied for 
flllSp . cf. |17j, by the following transformation of the unit square into itself 

x' = xy, y'=^^^—^. (116) 

1 — xy 

For example, the integral A°^\y^v3VA0 ^^^ ^^ calculated combining first the propagators D2 
and D^ by introducing a Feynman parameter, then performing the integral in k^. After this the 
result is combined with the remaining two propagators and the ki integral is carried out. One 
obtains 

<...3.40 = C f dx f dy f rf^x''+'^"*-3-i(l-x)-'^^+i+il/"+^+'^-i 
Jo Jo Jo 

X (1 - y)-^3+l+f^<^24-3-f (^ _ ^y-,-l 

X [z{l -X)+ X{1 - y)f-'''^'^+' . (117) 



Here we use the notation 



T^i,j...k = ^ ^i (118) 

l={i,j...k} 

and C = {-iy^'''i^r{ui234^ - 4 - 6){miy-''^^'^''+%A ■ p)"+^ In the case a = n,b = 0, after 
interchanging x and y, the integral flll7p is of the form given in flll4p . On the other hand, if 
a = 0, 6 = n, we obtain an expression of the form flllSp . The change of variables according to 
(ITT6|) yields 

<:..3.40 = C fdx' fdy' f dz x'-{x' + y' -x'yT^--'-'x'^^-'y'-^-^+'+'^ 
Jo Jo Jo 

X fl — x')~''i^^^**+'^+^^('l — ■y')^'^4+^+f ^^24-3-| 

^ {I - zr-' [y' + z - y' z]'-''''''^' , (119) 

which is of the form (11140 as desired. 

The easiest way to calculate the integrals with five propagators is to write them in terms 
of 4-propagator integrals, using integration by parts (IBP) identities [IE]. For the 5-propagator 
A-type integrals which appear in expressions (I106p to (I109p . one obtains 

An,0 _ 1 //.n,0 _ An,0 , .n,0 \ /.r,Q\ 

^11111 ~ 1,^12101 ^02111 ^ ^11102^ \^^^) 

AO,n _ \ ( /\0,n I 40,n i ^O.n _40,n _i_/lO,n _/.0,n \ H 9Tl 

^11111 ~ 1,^21011 ^ ^11120 ^ ^01121 ^10121 ^ ^11012 ^10112;' \^'^^) 
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.n,0 _ _1 

^21111 ~ 

.0,n _ _1 

^12111 ~ 



^22101 ^21102 ^ 1 _ ^ V ^12102 ^^ ^^03111 -^^13101^ 

AO,n /iO,n /iO,?i /iO,n 

~ ^12210 ~ ^02211 ~ ^12102 ~ ^22101 

-'- / ~ .0,n I r, /lO.n /iO,n , d AO,n AO,n 



;i22) 



T^ 1 _ I, -^^31101 T^ -^^30111 ^21210 ^ ^^20211 ^21102 

_AO,n _r,AO,n I O40,n _r,40,n _40,n \ (^9'^'] 

^21201 ^^11310 ' ^^10311 ^^01311 ^11202^ • K^^'^J 

For the ii^-type integrals one finds 

<iii = Y^ (-fe^iirif-^ + in + b)AliTu-' - EtLi + K2111 - <i02) . (124) 

In this way, the 5-propagator ii^-type integrals can be obtained from 4-propagator integrals of the 
same type, together with the previously calculated Aniii- I^ 01241) the factor of n multiplying 
the integral ^ii"u~ can be absorbed in the integrand using integration by parts 

n f dx x"'V(x) = f dx x"" [f{l)5{l -x)- fix)] , (125) 



where we have assumed that x"'f{x)\x=o = and f{x) is regular at x = 1, which is the case for 
the diagrams being considered. Finally, we have 

pa,6 _ 1 f,p.a+l,b-l pa,b _ pa,b T^a,b \ /.ofil 

-'^11111 ~ T , I "-'^11111 T-r]^2101 -'^02111 T^ -'^11102 i ■ {^"^^J 

This relation can be used recursively, that is, once -Fn'm is obtained from only 4-propagator 
integrals, we can use it to obtain -Ffnii and then -Fiim- More details on the way to obtain these 
equations can be found in Appendix IB.ll 

The results for all of the required 4- and 5-propagator integrals, appearing in expressions 
fll06H113"]) can be found in Appendix IB.5I All of these integrals were checked numerically for the 
first few moments using Tarcer [l9]. It turns out that all of the integrals can be expressed in 
terms of Nielsen integrals [50] , 

.l)n+p-l /•! ^^ 



SnA^) = ) n, , / -ln"-^(^)ln^(l-^^) (127) 

{n - ly.pl Jq z 

Un{x) = 5„_i,i(x) , (128) 

partly weighted by denominators 1/x, 1/(1 — x)^, 1/(1 + x)\ k,l < 3, as well as the distributions 
(5(1 — x) and 'Dk{x). 

The diagrams with two photons being attached to the operator vertex are more simple, since 
in this case only four propagators contribute. Furthermore, three factors of ZX occur and only 
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one term survives in the numerator after taking the trace. For diagram [S} we obtain 
• 1 



'5,j 



'o I ^ 



X 



- ln(l - x) - ln(a;) - 1 + (1 - C2)<^(1 - x) + 2)o(x) 

1 — X 



'Diix) 



X 



1 — X 



(42)o(a:) - 4)C2 - 4 - 62)2(x) - 42)i(x) + 42)o(x) 
(31n2(a;) + 21n(x) ln(l - x) -41n(x) - 2Li2(l - x)) 



+ 2 ln(l - x) + 6 ln2(l - x) + (10 - 4(2 - 6(3) 5(1 - x) 



-32 



+ 24C2 - 64C3 + 56C2 ln(l - x) + 48 ln(l - x) - 12 ln2(l - x) 

C2 ln(x) - 881,2(1 -x) + 5Li3(l - x) 



112, 3,^ , 8x 
^•^(l -x) + 



1 — X ^ 

- 3 ln(l - x)Li2(l - x) - 2 \n'^{x) ln(l - x) - 2Li2(l - x) 

- - In^(x) + 2 ln(x) ln(l - x) + 2 In^(x) - 2 ln(x)Li2(l - x) 



ln(x) ln^(l — x) 



32Li2(l - x) + (32 + 8C2 + 6i(3)1>o{x) 



112, 



- (56C2 + 32)2)i(x) + 162)2(x) + — D3(x) 



+ (152 - I2OC4 - 24C2 - 32C3)5(1 - x) } mf (A ■ p)" . 



:i29) 



iagrarr 


I [5] yields 








h, - 


Jo 


dx 


x"^ 


fl6 

16 

e 



- 1 - 2 ln(l - x) - X + 2!Do(x) + 2'Di{x) + (1 - C2) 5(1 - x) 
2x ln(l - x) + 3 ln2(l - x) - 2x - 2(2 + (2 + 2(2) 2)o(x) 



2Di(x) - 3D2(x) + (2 - 2C3) 5(1 - x) 



+ 32 - 64x - 64C3 



(20 - 12x)C2 + 16(1 - 2x) [ln2(l-x)-21n(l-x)] + 162)2 (a;) 

112 

— - ln3(l - x) + 56C2 ln(l - x) + (32 + 8C2 + 64(3) 2)o(x) 

/ 23 \ 112 

(32 + 56C2) 2)i(x) + ( 48 + 4C2 - -ttM 5(1 - x) + —V.^x) 



+ (-!)" 



16 / 2x \ 16 

— ln(x) + 1 - X ( 2x ln(x) + 2 - 2x 

e"' V 1 + X / e 



1 + x 



2C2 + ln^(a;) + 4 ln(x) ln(l + x) + 4Li2(-x) 
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64 - 64a; - 32(1 + x) \n{x) ln(l + x) + Ux \n{x) 



+ 



ix 



C2 \n{x) + 8 ln^(l + x) ln(x) + 8 ln(a;)Li2(-a;) 



l + x 
+ 4 In^(x) ln(l + x) + 8(2 ln(l + x) + 16 ln(l + a;)Li2(-a;) 

+ - ln^(a;) + 16Si,2(-a;) - 8Li3(-a;) - 8(3 ) - 16a; ln^(a;) 



^12 + 20x)C2 - 32(1 + a;)Li2(-x) 



The result for the unrenormahzed matrix element Aee is 



mf(A-p)". (130) 



i?>^ 



Si / dxx'^l — 



1 + 3x2 



ln(x) - 4(1 + x) ln(l - x) - 8 - 4x + 12©o(a;) 
4x2 + 2x - 1 



1 — X 



ln(x) — 7 + 3x 



1 — X 
+ 82)i(x) + (18-4(2) 5(1 -x) 

+ (l + x) M21n2(l-x)-iln2(x)-4Li2(l-x) -8C2) - 24T»2(x) 

+ (28 + 12x) ln(l - x) + 2 -— ln(x) ln(l - x) + (I6C2 - 2)2)o(x) 

1 — X 

- 40Di(x) - (IOC3 - 5C2 - 24)5(1 - x) 
1 + 3x2 



104 32 , , 

^-x + 16 + C2 

3 1 — X 



+ 



[6C2 ln(x) - 8 ln(x)Li2(l - x) - 4 ln2(a;) ln(l - 
1 + x2 



1 — X 
- (24-18x)ln(l -x) + 16 
64 



X 



1 — X 

51 



22 
+ ( — X + 32+ , ,^ . X, 

'3 3 1-x2 l-x 31-x3 



[2Li3(-x) - ln(x)Li2(-x)] 

16 \ , 9. . 80 



In^(x) + 



1 



3(1 -x) 



+ 56(1 + x)C2 ln(l - x) - (104 + 32x) 111^(1 - x) - -(1 + x) \n^{x) 

o 



178 



36x 



64 



140 



48 



3(1 -x)2 3(1 -x) l + x 



ln(x) 



+ 4^_!^_^ in(a.) in(i _ ^^ _ 2 ^ + ^^^ \n{x) \n\l - x) 



1 — X 



1 — X 



+ 32^^ [ln(x) ln(l + x) + Li2(-x)] - —(1 + x) ln^(l - x) 
1 — X 3 

5—117-2 "lO 

+ 4:^ i^ [ln(l - a;)Li2(l - x) - 11,3(1 - x) - 2(3] - 30x - — 

13x2 + 9^ ,^ , 4(16x2- lOx- 27) ^. ,^ 
- ^^^^^^^^^ - ^) + 3(1 -X) ^^-(^ - ^) 

+ (16 - 4OC2 + I28C3) 2)o(a;) + (24 - II2C2) 2)i(x) + 144D2(x) 
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224^ , , 
+ — 2)3(x) + 



67 
192 - — TT^ + (128 - 1441n(2)) C2 + 82(3 
45 



(5(1 -x) 



+ Ar 



24 
8 ) 5(1 - x) - 24'Do(x) 



\N 



4 / 1 + x^ 



(4 ln(x) ln(l + x) + 2(2 + 4Li2(-x) - In^(x)) 



e \ 1 + X 
+ 4x - 4 - 2(1 + x) ln(x) ) + 

4(x2 + lOx - 3) 



2(l-x)(45x2 + 74x + 45) 



+ 



3(l + x)2 
(C2 + 2Li2(-x) + 2 ln(x) ln(l + x)) 



3(1 + x) 

+ V^ { 8C2 ln(x) - 24C2 ln(l + x) + 36(3 - \ In'(x) 
1 + X \ 3 

- 4 \^{x) ln(l + x) - 24 ln(x) \^{\ + x) - 16Si,2(l - x) 



- 8 ln(x)Li2(l - x) - 48 ln(l + x)Li2(-x) + 40Li3(-x) 



-241n(x)Li2(-x) -48Si,2(-x) J + 4x 

2(9 + 12x + 30x2 _ 20a;3 - \^x^) 
^ ^,_ — — ln(x) 



;i + x)= 



3(l + x)3 
16 (x^ + 12x^ + 12x2 + 8x + 3) 



3(1 + x) 



Li2(l-x) 



-4x'' 



5 — X 



In'fx) 



mW^-vT 



(131) 



Notice the presence of the term proportional to A^, which is canceled by the counterterms. This 
term stems from diagrams with one-loop insertions, such as diagrams [Sji and [5)i. 

5.2 The OME A^^)'" 

The matrix elements AL are obtained from the fermionic one-loop insertions shown in Figure [7] 
supplemented by the corresponding self-energy diagrams. 





(a) (b) 

Figure 7: Feynman diagrams for the calculation of the massive two-loop operator matrix elements 

4(2),II 



One may use the fermionic one-loop vacuum polarization 

8e2 



nr(?) 



(47r)^/2 







{m\ — x(l — x)g2)2 -D/2 



(132) 
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in the corresponding one-loop diagrams, cf. Section |H The result for AL is 



4 (2), II 



tJ - I dOC X S ;r 



^(1 + x)-^2)o(x)- 3(5(1 -x) 



5x — 7 1 + x^ 



9 



3 1-x 



■ln(x) 



^(1 + x) ln(l - x) + ^Di(x) + ^2)o(x) - ^(2 + 3 ) (5(1 - x) 



76 572 

+ 27"""^ 
64 



4 8 

12x + - + 

3 1 — X 



ln(x) + 



128 



80 



32 



9(1 -x)2 27(1 -x) 
2 



9(1 -x)3 9 V(l-a;)^ (l-x)3 (1-x)^ 



ln(x) 



16 2('1 -I- t"^) 

+ -(1 + x) (ln(l - x) + ln^(l - x)) - ^^^_J \n\x) 

+ ^(1 + x)C2 - y (2)i(x) + D2(x)) I ml%A ■ pj^ \ . 



(133) 



5.3 The OME A^^)"! 



ee 



The flavor pure-singlet diagrams of Figure E] yield AeJ' ■ They can be calculated using the 
corresponding one-loop off-shell insertion of the diagrams in Figure 4, see |27] . 



crossed 





Figure 8: Feynman diagrams for the calculation of the massive two-loop operator matrix elements 

4(2),III 



For the insertion in diagrams [Hk,b one obtains 



4l^-^r(3 - D/2) r dx x"+^/2-2(i _ 3,)-2+D/2 

_2a;(l - x) .^r"^'7yi. (A ■ kf 



'^mlg^y 



:i-x) 



^-P + M2)3-0/2^ '^ (_A;2 + M2)3-^/2 

2(n -|- l)x — n k^/S.y + ky/S.^ 



{A-kf 



2 - D/2 (-P + M2)2-^/: 



+ (l-x 



n{l — 2x) — -Dx 



dfiv 



2 - D/2 (-A;2 + M2)2-^/2 



:iA-k) 
(A • kf 
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'1 — x)n- 



(n + l)(l-x)-l 



A^A, 



lSj'\k) 



1 - D/2){2 - D/2) (-P + M^y-^n 

, (A ■ kY'\^^ _ ^^2) r ^^ ^^ _ 3,2)D/2-l(^„ + (1 - x) 



(134) 



(47r)^/2 



X 



(A ■ k)/\^K - fc^A^A, 
(-P + M2)2-s/2 



(135) 



with M^ = m1/[x{l — x)] and k is the momentuiii of the photon hne. Using these results AL 
is given by 



4(2),III 



Jo 



1 — X 

3x 



(4x^ + 7x + 4) + 2(1 + x) ln(a;) 



4 
+ - 



5(l + x)ln2(x) - 



l + x 
3x 



^x^ - 17a; - 16) ln(x) + 



4(1 -x) 
9x 



(5x^ 



+ 23x + 14) 



2 1 

+ -(1 - x){Ax^ + 13x + 4)C2 + — (8x^ + 135x2 + 75x 

3x 



+ 32) ln2(x) + 



X 

304 80 



32 



x^ X + 108 



32 64(1 + 2x) 



9x 9 



l + x 3(1 + x) 



ln(x) 



224 1 — X 

x^ + 16 (x^ + 4x + 1) [2 ln(x) ln(l + x) - Li2(l - x) 

27 3x 



+2Li2(-x)] + (l + x) 



4(2 ln(x) + — In^(x) - 32 ln(x)Li2(-x) 
o 



161n(x)Li2(x) + 64Li3(-x) + 32Li3(x) + I6C3 
64 



182 



-x + 50 



32 800 



l + x 27x 3(1 + x) 



m''{A-py . 



(136) 



The factor [1 + (— 1)^]/2 occurs because of the different directions of the momentum flow in 
the external fermion lines in the diagrams of Figure El For all operator matrix elements, which 
are not the same in the unpolarized and polarized case factors of this type appear [51]. In the 
present case N is even. 



5.4 0{a^) Results 

From the pole-terms of the unrenormalized OMEs (1131(11331 11361) . resp. the renormalized OMEs 



4] 476|) . one may determine the 2-loop splitting functions Pe 



{i),i. 



x) to Pee (x) and the constant 
parts of the unrenormalized OMEs Tee ~ (x). Here the splitting functions result from a massive 
calculation. The operator matrix element ( ]13ip contains two branches. The contributions to 
the present process result from the first branch. The corresponding contribution to the 2-loop 
splitting function is found to be 

l + x^f 3ll 1 

<^ ln(x) ln(l - x) + - ln(x) ^ + -(1 + x) \n^{x) + -(3 + 7x) ln(x) 



5(1 

2^ 



X 



_(l_8C2 + 16C3)<5(l-x) . 



(137) 
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The pole part of the second branch, labeled by the factor (—1)^, is given by 



i^i)'^(x) ^ P£Hx) = ipW(x) 



1 + a; 



S2{x) + (1 + x) ln(x) + 2(1 - x) 



(138) 



where 



S2{x) = -2Li2(-x) - 2 ln(a;) ln(l + x) + - \n\x) - C2, 



(139) 



cf. go]. 

From the OMEs II and III the splitting functions 



P^^'^'^i 



-^Do(x) - ^ (1 + 8C2) 5{1 -x)- y 1^ + l{l + x) ln(x) - ^(1 - llx) 



l + x^ /I 5\ 2, 

-Inx + - + -(1 - 2;) 

1-a; V3 9/3^ ' 



(140) 



^l^'^^(^) - \pil'''\x) 



— I — X H — X \ mx 1 + xmxH 

2 2 3/ 2^ ' 9x 



(10 + X + 28x2) 



(141) 



are obtained in accordance with known results from Quantum Chromodynalmics [30] setting 
Ca = 0, Cp = Tp = 1. In the present calculation we choose Nf = 1. 

The unrenormalized operator matrix elements at the level of AL ( ITTHTSj) obey fermion 

number conservation, i.e. their first moment vanishes. The constant contributions to the un- 
renormalized OMEs (17TI473|) are given by 



f(i),i 



1 + 3x2 

1 — X 
'122 



+ 



[6C2 ln(x) - 8 ln(x)Li2(l - x) - 4 ln2(x) ln(l - x)] 

32 \ 1 + x^ 

X + 22 + ) C2 + 16^ [2Li3(-x) - ln(x)Li2(-x)] + 



1 — X 



80 



1 — X 



3(1 -x) 



'22 64 51 16 

+56(1 + .)C. ln(l -.)+(-. + 32 + 5i^—jj - ^— - 5i^—j3 



1 78 
-(92 + 20x) ln^(l - x) + ( 36x + 



64 



140 



48 



3(1 -x)2 3(1 -x) 1 + x 



In^fx) 



ln(x) 



1/. M 3/ N .a;2-8x-6, , ,, ,^ , ^1 + ITx^ , ,, 2 

a + x)ln^(x)+4 1„^™M„n ^\ o 1„/™M„2 



3 
112 



1 — X 



■ln(x)ln(l -x) -2- 



1 — X 



ln(x) In (1 — x) 



1 + x 62 
,1 + x) ln^(l - x) + 32^^ [ln(x) ln(l + x) + Li2(-x)] - 22x 

O 1 X o 

-4 . ^ Si,2(l - x) + 4-^ [ln(l - x)Li2(l - x) - Li3(l - x) - 2(3] 



1 — X 

4(16x2 -lOx- 27)^^ 
3(1 -x) 



1 — X 
Li2(l - x) + 14(x - 2) ln(l - x) 



224, 



+ (16 - 52C2 + I28C3) 2)o(x) + (8 - II2C2) Di(x) + 1201)2(0;) + —'Dsix) 
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+ 

+1 



^-|-^+(y-481n(2))C. + 58C3 



S{l-x) 



( 2(1 - x)(45x2 + 74x + 45) 2(9 + 12x + SOx^ - 20x3 - 15x 



+ 



3(l + x)2 
4{x^ + lOx - 3) 



3(1 + x) 



+ 



1 +x- 
1 + x 



2 r 



3(1 + x) 

(C2 + 2Li2(-x) + 2 ln(x) ln(l + x)) 

2 



ln(x) 



8C2 ln(x) - 24C2 ln(l + x) + 36C3 - 5 ln^(a;) + 40Li3(-x) 



-41n2(x)ln(l + x) - 24 ln(x) 111^(1 + x) - 241n(x)Li2(-x) 



- 48 ln(l + x)Li2(-x) - 8 ln(x)Li2(l - x) - 16Si,2(l - x) 
16(x^ + 12x3 + 12x2 + 8x + 3) . 



-48Si,2(-x) 



3(1 + x) 



-Lio(l-x) 



+4x r:^ In^(x) 

(1+X)3 ^ ' 



(142) 



f(i),n 



76 572 

27^ ~ ^ 

32 / 1 



12x + - + - 
3 1 



-^ 



X J 



ln(x) + 



128 



1 -X 



'l-x) 



+ 



'l-x) 



9(1 -x)2 
ln(x) 



+ 



80 



64 



27(1 -x) 9(1 -x)3 



2(l + x2), ^, . 4 



3(l-x)^^(^) + 3^^ + ^)^^ 



H (1 + x 

3 ^ 



(ln(l - x) + 111^(1 -x)) + (^- ^(2) Do( 

32 /S 1 41 1 \ 

-- (Di(x) + V,{x)) + ( -C3 - IOC2 + ^^j 5(1 - x) 



X 



(143) 



fW'™ = -(l-x)(4x2 + 13x + 4)C2 + — (8x3 + 135x2 + 75x + 32)ln2(x) + 50 
X 3x 



+ 



304 80 2 32 _^ 32 64 1 + 2x 

x^ X + 108 r -r^ 

9x 9 3 1 + x 3(l + x)3 



224 2 182 

m X X X 

^ ^ 27 3 



+ 16i^(x=^ + 4x + 1) [2 ln(x) ln(l + x) - Li2(l - x) + 2Li2(-x)] 



3x 

+ (l+x) 



800 

27x 



14 , 
4C2ln(x) + —W{x) - 321n(x)Li2(-x) - 161n(x)Li2(x) 



+64Li3(-x) + 32Li3(x) + I6C3 



32 



+ 



64 



1 + x 3(1 + x)^ 



(144) 



The constant contributions (I142H144I) enter the single differential cross section hnearly and all 
contain terms of the kind ~ 1/(1 — x)^ and ln(x)/(l — x)^ etc. The massless Wilson coefficients 
of the Drell-Yan process to 2-loops are free of such terms, cf. [201121], like all the other Mellin 
convolutions which contribute. For the processes II and III non of these terms were found in pQ 
resp. [25], while terms like In (x)/(l — x)^ are present in the cross section for process I. We 
therefore conclude that the use of standard matrix elements for local operators alone, f lT71 - [2T|) . 
are not sufficient to reproduce the constant terms. This aspect requires further investigation. 
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It occurs for massive external fermion lines in contrast to the case of massless external parton 
lines. There the constant terms are correctly reproduced in the limit s ^ m'^. 

Finally we present the O(aQ) contributions up to the logarithmic term L in ( 140] - [42|) . which 
have been computed based on the present 2-loop calculation of the massive operator matrix 
elements with an external massive fermion line. They are given by 



T, 



2,Ia 



16an 



4Di(x) + 3Do(x) + {l- 2C2I 6(1 -x)- 2-''^"'^ 



1 — X 



(1 + x) ( 2 ln(l - x) - -(ln(x) - 1) ) - 1 + x 



-8Di(x) + (-7 + 4C2) 'Doix) - ^(11 + 7x) ln(x) 



1 _)_ 3;2 / J^]^ 

j — In (x) H ln(x) + ln(l — x) ln(x^ 



+ (1 + x) ( -2C2 + 41n(l -x) + - In'(x) ) + 3 + 4x 



+ |-^ + yC2 + 3C3)<5(l-x) 



(145) 



T; 



2,IIa 



IQan 



iDo(x) - 1(1 + x) + 1^(1 -x) 



+ 



1 +x^ 

1 — X 



- ln(l — x) ln(x) 



9 



:i-x) 



17 
12 



5(1 -x) 



L } (146) 



T- 



2,IIIa 
ee 



16an 



1 1 1 

3x ^ 4 ^ 4' 



1 + x) ln(x) H 1 X x^ 



+ 



;i + x) (2 ln(x) ln(l - x) - ln2(x) + Li2(l - x)) 



4 4 \ / 2 14 

"^ ' 3x ^ ^ ~ ^ ~ 3^ / ^^^^ ~ '^'' ^ V 3x ^ ^ ~ 2^ ~ 3^^ ' ^^^^^ 



\ — X H — X 

9x 3 3 9 



(147) 



These terms agree with the corresponding terms in Ref. [T], Eqs. (2.30, 2.40, 2.42, 2.43) and for 
(IT16|) with Ref. [25]. 

6 Conclusions 



We have calculated the QED matrix elements with an external massive fermion for the local 
operators of the light-cone expansion to 0{a'^) up to the constant term in the dimensional 
parameter e in the MS scheme. The OMEs can be renormalized applying the technique having 
been developed recently in [27]. The renormalized OMEs obey fermion number conservation. 
Various technical details of the calculation and intermediate results are provided. We investigated 
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the factorization of the Ola"^) QED initial state corrections to e+e" annihilation into a virtual 
boson for large cms energies s ^ ml into massive OMEs and the massless Wilson coefficients 
of the Drell-Yan process adapting the color coefficients to the case of QED, as being proposed 
in Ref. [1]. We have shown by an explicit calculation, that the representation works at one- loop 
order and at two-loop order including all terms up to the linear order in ln{s/ml). In the constant 
term in 0(0^) the OMEs bear a few structural terms, which have not been obtained in previous 
direct calculations [ll|25], despite a large part of other terms appear as being expected by the 
factorization. Further studies are needed to reveal the reason for this. This finding appears in 
contrast to the case of massless external fermion and boson lines, where the corresponding cross 
sections have been shown to factorize including the constant terms in 
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A Convolutions of one— loop quantities 



In Eqs. fl40] - |42|) a series of convolutions of one-loop functions occurs. Some of them where given 
in the Tables in [T5l|52]. |j A few more convolutions are needed which are given by : 



2)i(x)®2)i(x) 

Di(x) (g)ln(l -x) 
Di(a;) (8)a;ln(l - x) 



V,ix) - 2C2Di(x) + 2C32?o(x) - ^6(1 - x) 

2 [Li3(x) - Cs] - ln(x)Li2(x) - ln(x) ln2(l - x) 
1 — X 
Li3(l -x) + 2Li3(x) + [ln(l - x) - ln(x)] Li2(x) 

+ ^ln=^(l-x)-2C2ln(l-x)-C3 

xLi3(l — x) + 2xLi3(x) — [1 + x [ln(a;) — ln(l — x)]] Li2(x) 
+- ln3(l -x) + {l-x) ln2(l -x) -{1-x) ln(l - x) 
— a;ln(a;) — (1 — x) ln(x) ln(l — x) — 2x^2 ln(l — x) + x(C2 
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g>x^ 




ln(x)x ( 


g>x^ 




ln(x) C 


g>X^ 




ln(x) 
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(8)x 




ln(x) 
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® 1 


ln(x) 


®ln(l 
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ln(x) (g 


)xln(l 


-x) 


xln(x) 


®ln(l 


-x) 



^[31n(x) + l-x3] 



■ x\ 



X [ln(x) + 1 

- [2 ln(x) + 1 - x^] 

^[21n(x) + l-x^] 



■ X 



- [ln(x) + 1 

X 

Li3(x) -C2ln(x) -Cs 

- [Li2(x) - C2] + (1 - x) [ln(l - x) - 2] - ln(x) 

ri 



X ln(x) X ln(l — x) 

ln(l — x) ® ln(l — x) 

ln(l — x) ® xln(l — x) 

xln(l — x) ® xln(l — x) 



X [Li2(x) — C2] ^ (1 ~ 3;) ln(l — x) + X ln(x 
[Lis (a;) -C3-C2ln(x)] 



-ln(x)-l 



- {25i2(x) + 25i2(l - x) + ln(x) [Li2(l - x) - C2] - 2(3} 
'' ' ^^'- '' x) + ln^(l-x) - "-'" 

[Li2(x) + ln(x)] 



(1 - x) [Li2(l - x) + ln^(l - x) - ln(l - x)] 

+C2(2x - 1) - X [Li2(x) + ln(x)] 
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148) 



(149) 



-Cs) 
(150) 

(151) 

(152) 

(153) 

(154) 

(155) 

(156) 
(157) 

(158) 

(159) 
(160) 

(161) 

(162) 



We finally list some Mellin convolutions of splitting- and related functions : 



\ [Pf^ ® Pf^) (x) = 642)i(x) + 48Do(x) + (18 - 82(2) 5(1 



— X 



ln(x) 
-32^^^- 16(1 + x) 



16 



1 — X 

l + x^ 

1 — X 



3 3 

21n(l-x)--ln(x) + - 



16(1 -x) 
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(163) 



^The distribution- valued contributions to ()148p were also given in [5 3) 
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UP^>Pi^)(^) = 16 
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'l + x)ln^(x) + 8xln^(x) 
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-8 + 8x x" 

9 X 9 



64 1 



In f 1 — x) 
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+ + + X X" (172) 

27x 9 9 27 ^ ' 

^2)3(x) + 722)2(x) + (24 - 48C2)2?i(x) + (64(3 - 2OC2 + 32) Do(x) 
+4(1 + x) I 4Si,2(x) + 6C2 ln(l - x) - I2C3 ln3(l - x) - 3 ln(l - x] 

-2- C2 ln(x) - 16^ ln(x) [ln(l - x) + ln^(l - x 



1 — a; 1 — X 
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B Integrals 

In the following we present some decompositions of five-propagator integrals in terms of four- 
propagator integrals using integration by parts [Hj. In some cases Mellin-Barnes representations 
are applied for checks. In a sample calculation we illustrate main steps of the computation of 
the massive OMEs. Furthermore, we list a series of complicated integrals and finally summarize 
the structure of the individual diagrams in terms of Nielsen integrals. 

B.l Integration by parts 

Eqs. f ll20H123"|) are obtained from the integration by parts relation 






d^h d'^k^ d f (A ■ A;i)"(A ■ k2Y 

where q^ = k^, kt^ and l^ = p'^, /cf , /ca • Fo^ example, (11201) results subtracting 8.i{iii{k2, ki) from 
Siiiii(^2; ^2)- On the other hand, the most complicated equation of this family, namely, fll23p 
can be obtained from 

^21111(^15 ^2) — £21111(^1? ^1) + Sii21l(^l5 ^2) — £11211(^1? ^1) 

-(D-5)[£?2';n(fci,A:2)-e?2m(fci,^i)] = 0. (175) 

In the same way, it is possible to obtain relations for v4"{j^^]^ and A^^^-^-^^, which are needed to 
obtain the i?-type integrals from (1124p . 

Sometimes it is even possible to use (I174p to reduce 4-propagator A-type integrals to a linear 
combination of integrals with just three propagators. For example, it can be shown that 



An,0 _ 1 j ATifl _ ^ 
^31110 ~ 1 1 ^ ] ^32100 



^23100 ^ 1 _ ^ 1,^14100 ^04110;' 



(176) 



It turns out to be much easier to calculate ^aino from this equation than by introducing Feynman 
parameters directly. 

Equations for the 5-type integrals can also be obtained by setting one of the powers of the 
propagators in (I174p to zero. For example, from 

S?mo(^i.P) + £?mo(^2,p) = , (177) 

one may obtain 

-D12IIO — 2^ Pj^lUlO ^ ^21110 r, 01210 + 2 21010 ' {^'°) 

and from S.i{iQi{ki, k2) = 0, one may get 



Tjn,0 _ 1 r ^n-1,1 _ Aufi , An,0 _ An,0 , An,0 _ /. 

-^21101 ~ 2 L 11101 ^21001 ^ ^21100 ^20101 ^ ^11200 ^ 



n,0 
10201 



~^iioo2 + ^10102 + 2"^ ^11102] • (179) 

These kind of equations can be used to check the results obtained directly using Feynman 
parameters for this type of integrals. 
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Let us consider now the ii^-type integrals. Using 

n-l 

A-k2- A-ki 



£(A . fc.v(A . k^r-^-' = '^;,':i":'f;/-'= 



3=0 

we can rewrite them as 



rpa,b ja,n+b jn+a,b 



where 



j: 



,, , d^k, d^k2 iA-k,riA-k2) 



uiu2u,u,us I (27r)«(27r)«Z}A^r^2'^3'^r^r ' 



'5 

with D^ = A ■ k2 — A ■ ki. We define also 

r^a,b ^ f d^'k, d^k2 (A ■ A:i)"(A ■ k^)' 



Then, from 



d^d^^ f _ {A-k.riA-k^Y 

{27r)D {27t)d dkl^ V ' DaD,D2D,D,D, 



and using 

and 

Eq. f ll24p is obtained. Similar manipulations lead to (IE 



T^a+l,b T/-a,b+l ja,b 



ja+l,b ja,b+l A^fi 

'^ V\U2VZV4,Vz '^ UlU2US,ll4,Ur, U\l'2V3U4,V^ 1 



180) 

181) 
182) 

183) 

184) 

185) 
186) 



B.2 Mellin-Barnes representation 

The 5-propagator integrals can be also calculated and/or checked using a Mellin-Barnes repre- 
sentation [51]. To do this, we follow the ideas developed in [5S], where a Mellin-Barnes repre- 
sentation was used to study the massless two-loop two-point function with arbitrary powers of 
the propagators. In fact, these ideas have been already applied to the calculation of massive 
operator matrix elements with external gluon lines [56] . 

The presence of the mass and the Mellin variable n make our integrals more complicated 
compared with the massless two-point function, so in order to obtain the Mellin-Barnes rep- 
resentation of lowest dimensionality, it is necessary to choose the right momenta fiow in the 
diagrams. Let us consider the integrals At;^y^yo,UAV5- ^^ turns out that changing variables 

ki^k2 + p, k2^ki, (187) 

so that now the propagators in (IIOOI) become 

Di = {k2 + pf -n? , D2 = kl-m\ D^ = ^ 



D4 = {k2 -ki+ pf , D^ = {ki - k2 



2 ! 

2 _2 



one can obtain a two-dimensional Mellin-Barnes representation. This is achieved combining 
the propagators D2 and D^ with a Feynman parameter, and then combining the result with 
D4 introducing another Feynman parameter. After completing squares and performing the ki 
integral, one obtains 



' '^ <■ fk^ . r(z/245 -2-5/2) (-l)'-^^^ 
v^,v2,v,^v4.v, I —I -. I ^2;r)^ ' r(z/2)r(i/4)r(i/5) (47r)^/2 



^^^L....u4... = 1^ dxj^ dy I :f,^^ 



31 



^3 



D'i'Dl' [yml - x{l - y)Di- xy{l - x)D:i^''^^'-^'^ 



Now, the denominator can be split twice using 



1 Y- 



/+joo 
da X"F-"r(-CT)r(a + z/) , (189) 



(X + YY 2mT{u) 
and we are left with only two propagators, which can be calculated directly. One obtains 

^uuu2,us,u,,u, ^ r(^,2)r(i/4)r(i/5) ^^ ' {h-k)\k\{2mY 

X r^. r-~" r(-.)r(-.)£<^±i±i^)I<^±l±^ 

J-ioo J^ioo r(20- + T + k + V2b) 

/ X r(-(T - r + z/13 - 2 - 5/2) 

r(z/3 - cT)r(z/i 



X 



r(2o- + r + A; - z/i - 2z/3 + 4 + £:)r(-cr + a + 6 - A; + z/3) 
r(a + r + a + 6-z/i3 + 4 + £) 
^ r(r + & - A; + z/4)r(-r - i/as - 2^/4 + fc + 4 + e) 

r(6 - z/245 + 4 + £) ■ ^ ^ 

For this expression, we can use the package MB [57] to check the integrals numerically, which can 
be done for up to relatively large values of n. This is the main advantage of this method over 
Tarcer [l9], which allows to check only the first three or four momentq^. One may try also to 
obtain analytic results starting from equation ( I190p . but unfortunately, unlike the case where we 
have an external massless particle [56], this turns out to be rather complicated, although being 
possible in some cases. 

This method can also be applied to check the 4-propagator integrals. For example, the 
denominators in flll7p and f lllQp can be split using equation fll89p just once, leading to a simple 
one-dimensional Mellin-Barnes representation. In spite of this simplicity, on occasions the MB 
package cannot find a proper contour. This usually happens for integrals with high powers of 
the propagators, like the 4-propagator integrals appearing in equations f ll22p and fll23p . which 
have a highly singular structure. As it is well-known, this problem can be cured introducing an 
additional regularization parameter. 

B.3 Sample calculations 

As it was mentioned before, the method of integration by parts has the advantage that it allows 
to write all 5-propagator integrals in terms of 4-propagator ones, which are simpler to calcu- 
late because they can be parameterized using only three Feynman parameters. In the case of 
Eq. f ll23p . this simplicity is somewhat spoiled by the high powers of the propagators in the 
4-propagator integrals. Let us consider one of the integrals appearing in that equation, namely. 



^?mo ^ f dx f dy f dz x'^^^il - xY'^^^'ix + y - xyy-' 
Jo Jo Jo 



^''Tarcer can be modified to handle larger values of n on the expense of very long computational times. 
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X y-^+''^{l - yy/h-^-'/\z + y- zy)-^+' . (191) 

It is not difficult to see that the singularity structure in e is shared by the three Feynman 
parameters in a way that is not easy to disentangle. If one tries to expand in e directly, the 
integrals will be ill-defined, and it is difficult to find suitable subtraction integrals to cure this. 
To solve the problem, we decomposed the integral as 

A?f3ioCx/i + /2-/3, (192) 



where 



/i = -- [ dx [ dy x"+2(l - x)-^+^%x + y- xyf-' 



e 



^0 



\l-e 



xy-^^^'\l-yr'\ (193) 

h = f dx f dy f dzx'^+^il-xY^+^'ix + y-xyf 
Jo Jo Jo 

X i/-^+^/2(l - yy+'/^z-'/^{z + y- zy)-'+' , (194) 

dx f dy f dzx^-^^il-xy^+^'ix + y-xyY-' 
Jo Jo 

X y~^+'/\l - yy+'/^z-'/\z + y- zy)-^+' . (195) 

Here, integral Ji was obtained performing integration by parts in z. It turns out that the change 
of variables (I116p . which was used before to write the integrals as Mellin transforms, can be used 
now to obtain expressions that can be safely expanded in e. For example, for the integrals I2 
and I3, we change variables by 

y = y'^', z = y}^ 'f) , (196) 

1 — z'y' 



which leads to 



I2 = I dx f dy' f dz' x"+=^-^(l - x)-=^+2. A ^ L^y'/ 



Jo JO 



X 



X y'-^+'z'-^+'^'{l - z'y^'il - y'z'Y , (197) 

/3 = f dx I dy' I dz' x"+3-^(i _ a;)"3+2. A ^ \jl^y'z'\ 
Jo Jo Jo \ ^ J 

X y'-^+'z'-^+'/^l - z')~'/\l - y'z'Y . (198) 

Now, the expansion in e (not including the term (1 — x)^^^'^^) will produce only logarithms that 
are regular at x = 1, and can be calculated using 

rl -1 r rl „ -a+l+be/1 _ Ace 

/ dy y~^^'^il ~yr In' (l + xy) = z — - -kx dy^ — ^ li^ln'-\l + xy) 

Jo 1-a + be I Jq I + XV 



+ce [ dyy-'^+'^''il-y)-'+^'\n'il + xy) 
Jo 



(199) 
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recursively to analytically continue the expressions by reducing the highly singular powers in the 
integration variables. A few other integrals which appear in (11231) can be performed in the same 
way. 

The E- and F-type integrals can be written in terms of Feynman parameters from the 
corresponding expressions for the A-tjpe integrals using 

n-l 
rpa,b _ V^ Aa+j,b+n-l-j (OQQ) 

j=0 
n-l 
Tra,b _ \^( A . j^)ri-l- j Aa+j,b (2Q^^ 

U1U2U3U4U5 Z-^^ ^' T^lJJ2VZT^A^5 ' K'^^ ^ J 

For example, 

„i »i „i ii,—± 

Jo Jo Jo „_n 



j=0 



^l n—1 



1 .1 .1 

dx dy dz 
Jo Jo 



j=0 

X z^^-'^'^[yz + x{l - x)(l - y^Y 

^n^i-en y-^l\l-y)-^l'z-^-^l' 

{x + y — xy)[xz + y{l — y){l — x)^]^^ 
n{x + y- xy)-^-'{l - xfy^-'/^{l - yY^l^z'^-^l 



[xP' + y — x'^y)[yz + x^(l — y)y 



(202) 



and 



where 



and 



„1 „i „i n-i 

Ekuo oc dx dy d^5^x"+i-^-^/2(l-x)^/V/'(l-l/r+' 

Jo Jo Jo j^Q 

X z-^"'/^il - z)[z{l -x)+ xy]-^^' 
= la-h, (203) 



"1 /•! /•! -l+e/2(i_\2-e/2-l-e/2(i_\ 

h= I dx I dy I dz x-^^ y JL1\ n\.n ^ Jl ^ (204) 



JO JO 



{l-x)-^-/^[zy + {l-y){l-x)Y 



,1 /•! /•! y£/2|'i _ WI+e/2 -l-£/2|'i _ N 

h= dx dy dz x"+i^ ^,/, , ^^ , , ^ \, . (205) 

Jo Jo Jo [l-xY ^''[x + y-xyY[y + z-yzy '^ 

Integral la is particularly difficult. One way to perform it is to use the integral representation 
of the hypergeometric function |58J 

^F,{a,p-r,z)= ^ / dtt^-\l-tr<-^-\l-tz)-^ , (206) 

^[P,l - P) Jo 

for Re 7 > Re /3 > and then use the following analytic continuation [59] 



2Fi(«,/3;7;z) = (l-^)-^2i^i(/3,7-a;7;^--) , (207) 
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to write 
with 



hi 

Ia2 
Ia3 



la — lal + Ia'2 ~ hs ) 

l(l-x)-l+=2/"l+^/V/2(l-z)- 



/•l /•! nx 

K dx dy dz x^' . , . . ,^ 

h k Jo [1 - X + {x - z)y]--/'^ 



(208) 

(209) 
(210) 
(211) 



where K = —- + eC,2 + yCs + ^'(e^). The integral lai, can be performed straightforwardly. For 
integral hi we change variables z = xz', while for 7^2 one substitutes z = {1 — x)z' + a;, which 
leads to 



Jo Jo Jo 



[1 - X + (1 - z')xy]-'/^ 



'a2 



1 /•! 



K dx dy dz' x 



J rr.n+l 



'l~xY'^y-^+''\l 



z 



JO JO [{l-x)z' + x\--l^[l-z'y]--/'^- 



(212) 
(213) 



Now integral Ia2 is easy to obtain, and integral hi can be done using another analytic continu- 
ation of the hypergeometric function, namely 

r(7)r(/3-a) / iv 



.F,(a,/3;,;.) ^ ™^ .F, («, « + 1 - ,; a + 1 - /3; i 

, r(7)r(a - /3) / ly / 1 

+ r(a)r(7-/3) l"J 2F,^/3,/3 + l-7;/3 + l-a;- 



(214) 



which yields 



lai = K I dx j dy I dz 
Jo Jo Jo 



X (1 - x)-^+'y-^-'z'/^{l - z] 



■n+2+s/ 



e/2 



X (l-xz)'' 1 + 



I — X 

(1 — z)x 



e/2 



+ {l+lc2+jC3]x-^\i-xr^'^ 



.-/2 



xz''^ l-z)-'^' 1-xz 



e/2 1 



(215) 



This integral can now be computed using known integrals and specific integrals given in the 
following Section. 

B.4 Polylogarithmic integrals and analytic continuations 

In the following we list a series of integrals over polylogarithms derived in the present calculation 
beyond those which were given in pOj. They may be of use in other higher order calculations. 

"^y T^JT^^'' i"^) = 2C2ln(x) + ^ln3(a:)+41n(x)Li2(-x) 

+3 ln(x)Li2(x) + 2Li3(x) - 2Li3(x2) , (216) 
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X 



"^^ (1 - (1 - ^)y) 



-Lio 



i-y 

x'^y 



X 



dy ^ 
l-xy 



hi-} 



1 — X 



x'^vO'-y) 



X 



^ \{l-y + xy){x + y-xy) 



-Li2(l -x) - -In^ix) -C2 , 

2Si,2(x) + 4Si,2(l-a;)-4C3 

— \n(x) ln^(l — x) , 

1 + x 



(217) 



(218) 



1 — a; 
-41n(l - X 



[Li2(x2) + 21n(x)ln(l + 
2x 



a; 



I — X 



[C2 + 21n(x)] , (219) 



T* 1 

dy ln(l — y) ln(l — xy) = Li3(l — x) + - ln(a;) ln^(l — x) 

1 — xy 2 

-^ln=^(l-x)-C2ln(l-x)-C3 



x{l — x) 

(l-a;i/)2 
1 



ln(l — y) ln(l — xy) = Li2(a;) + ln^(l — x) + ln(l — x) 



dy z 

l-xy 



ln{y) ln(l - xy) = Si,2(x) , 



T* 1 

'^1/7^; 7^1n(l/)ln(l -xy) = Li2(x) + -ln^(l -x) + ln(l -x) . 

[l-xyy 2 

In all of the results given above, it is assumed that < x < 1. 

Some double integrals that were also used are shown below, where m = (1 — x)/x. 



(220) 
(221) 

(222) 

(223) 



^ [^ zlnil — zy) In (1 + uyz) 1 + u 

dy / dz ■_ = 2 

Jo 



1- zy 



u 



1 



ln^(l + m) - ln(l + m) + Li2(-M) 



2 + 



ln(x) — Li2(l — x) 



1 — X 



(224) 



\ /"^ zlnil- z)\n{l + uyz) l + u 

dy I dz = 3 + 

Jo 



1 - zy 



u 



- \n\l + u)-2 ln(l + u) + U2{-u) 



3 + 



21n(x) -Li2(l -x) 



1 — X 



(225) 



Jo 



1 



1 - zy 



1 



4-2^ + 2^-^'^^^^ + ^ 



----A )ln(l + M) 



l-2x 
2(1 -x)' 



In'(x) 



1 -4x 
2(1 -x): 



ln(x) 



7x-l 



4(1 -X 



(226) 



\ /"\ y2;^ln(2;) ln(l + -uy2;) 
dy I dz 

Jo 



I- zy 



-^ + \-'^S,^^(-^) + U,{-u) 



+ 
+ 



1 1 , . 
H ln(l + M 



2%^ u 
■7 3 _ 1 



Li2(-M) 
ln(l + u) 



5x 



5 



4(1 -x) 2 
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-C3-ln(x)Li2(x) 



— In^(x) ln(l — x) + Li3(a;) — Li3(l — x) 



+ 



X 



3x^ 



1-x 4(l-a;)2 4 



+ I-C2] ln(a:) 






Jo 



1- zy 



3 - ln(l + u)U2{-u) + ln(l + m)C2 + -U2{-u) 



-2(1 + -) ln(l + u) + Lisi-u) - 2Si,2{-u) 



3 + 



1 — X 



1 X 



-C2 ) ln(x) - 2^3^^"^ *^^ 



X 



1 — X 



io 



1 - zy 



Li2(l - x) - Li3(l -x)- Si,2{l - x) (228) 



1 



u 



1 + 1 + - - C2 ln(l +u) + ln(l + u)U2{-u) 



-Li3(--u) + 2S'i,2(-M 
1 



-1 + C: 



ln(a;) 



1 f\ zln'jl+uyz) 

ay / az ; 

Jo 1-^1/ 



I — x^ 
+Li3(l-a;) + ^1,2(1- a;) 

2+ ( 1 + - j [\n^{l + u) -21n(l 

ln(x) 



(229) 



u 



2 + 



\n^(x) 



+ 2- 



1 — X 1 — X 



(230) 



^7 /"^ 7 y-2^1n (2;) In (1 — xy) 
ay / dz — 

Jo 



1- zy 



^^f^^^zHn{l-xy)\n{y) 
Jo 



I- zy 



',,/"<;, £!if?(l-^') '"(!-!') 

Jo 



I- zy 



15 1 / 1 2 V ^ 1 V 

— a; \ X X Co H — x G 

94 2V 2 y^^2 ^^ 

-ix^ [Li3(x) + 5i,2(a:)] - ^ (1 - a:') ln2(l - a:) 

+i(l-x2)ln(l-x)C2 (231) 

5 4,13^ A 2 1 a;3\ 

X X ^ — X Cs + \ -X H — X C2 

24 9 3 ^ V3 6 9 y ^ 

+^(l-x=')ln(l-x)[Li2(x)-C2] 

y-^ + 01n(l-a:) + i(2-x^)5,,(x) 

-l(l-2x=^)Li3(x) + l(l-x3)ln^(l 

— x^ H \ Li2(x) 

3 V 2 3; ^ ^ 

I + |(1 + x)C2 + -x\z - -(1 - x) ln(l - x) 

-- (1 - 2x3) Li3(a;) + - (1 - x^) ln(l - x)Li2(x) 
o o 
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X 



(232) 



Jo ^-zy 



Jo ^-zy 



Jo ^-zy 



1 .1 

dy / dz 
Jo 



z^ In (1 — xy) 
I- zy 



\ /"\ z^ In (1 — zy) In (1 — xy) 

dy I dz 

Jo ^- zy 



+ M - '^xA Si,2{x) + ^(1 - x)U2{x) 

+ ^{l - x) ln2(l -x) + -{l- x^) ln^(l - x) (233) 
12 9 

/ 7 a; 1 \ , /-, X fx x'^\ . 



+ fi^ - i) U^(x) + (1 - f ") Su(i) 



X 



+ - (1 - x^) ln(l - a;)Li2(x) (234) 

2 
-o + ^^C2 + (1 - a;^) Li2(a;) ln(l - x) 

+ - (l--j \w(\ - x) ^ Sx,2[x) ^ x^\^\^{x) 

(\ \ x^ 

+x(x - l)Li2(a;) +7-3; + — + — 1 In^fl - x) 
V 3x 2 



+ l(l-x^)ln^(l- 



a; 



(235) 



■H 2 



— X 



^ iXj I v/y / iXj ti-J ti-J \ , tXj . 

T2+^+U"^"6j^' + T^' 

+ ^ (1 - x'') ln(l - x)C2 - y [5i,2(:r) + Li3(x)] 

-i(l-x^)Li2(x)-l(l-x^)ln^(l 

— + — a; + -a;^ ) ln(l-a;) (236) 

12 36 9 y ^ ' ^ ' 

x^ 2 

y + a:^C2 + 3 (l - ^^) 1^(1 - x)\^i2{x) 

--(1 - x) ln(l - x) + - (1 - a;3) ln^(l 
o y 

+ ( -^ - -a;2 + x^ ) Li2(x) + -x^Li3(x) + -5i,2(a;) 



X 



3 3 

x^ + -x^ I ln^(l - x) 

2 3 3 2;^ ' 

3x 11x2 / 3 5x llx2\ , ^ 

— + ^7^ + I -T + — + —TT^ 1 ln(l - x) 



(237) 



+ 



18 ■ V 4 36 
1 X x^ llx'^ 

\ 

3 



9 6 



, Li2(x) H x^C2 

18 / ^^ ^ 18 ^ 



y(l-x=')[Li3(x)-5i,2(:r)] + yC3 



+ (1 - x^) 
38 



— ln3(l-x) + -ln(l-x)Li2(x) 



I X X llx" 

^ ~18 ~ 12 "y^ ^6" 



In^d -x) 



(238) 



^ r^^ yz^ ln{l - zy)ln{l - xy) 
Jo ^-zy 



7 llx 



9 12 



a; — ^r^ — ^^ 1 ln(l — x) 



12 



11 n 3 a; 
+ I —x^ - - 

24 8 



36 9a; 



1 



^ +— lln^fl-x) 
4 ox 



io 



1-2:1/ 



+— a;V2 + — (1 - x^) In^fl - a;) 
12 12 ^ ; V ; 

+ \ x^ Li2(a;) + -ajYs 

VSa; 4 2 12 ) ^^ ' 2 ^^ 

+ ]^{l-x^)[S,,^{x)-U:,{x) 

+-(l-a;2)[ +ln(l-x)Li2(a;)] 
5 a; /_, a;\ a;^ 

r-2i^-2^^+Y^^ 

-- (1 - x) ln(l - a;) - i (1 - x^) U2{x) 
+\ (1 - x^) ln(l - a:)C2 - ^ (l - x^) \n\l - 



(239) 



X) 



1 



a;^ [Li3(a;) + 5*1,2 (a;)] 



(240) 



io 1 - ^y 



3 1 

'-X- - (l-a;^)ln(l - a;)C2 



+ \{l-x)Hl-x)-^-(l + ^-)u,{x) 

X \ ^ , -. I 2 Iit; /„\ I 2, 



+ ( 1 - y 1 '5i,2(a:) + ( ^ " 3 ) ^13(2^) + 3^ Cs 



i (1 - a;^') ln(l - x)Li2(x) 



(241) 



^^ r^,^£!Ln(l-,)ln(l-a:,) 
Jo 



I- zy 



a;C2 + 3:^3 + ( 2 ~ ^^ ) '^i'2(a;) 

+ (1 - a;2) ln(l - a;)Li2(a;) - ( - - a;^ J Li 



^{l-x^)\n\l-x) 
6 



(242) 



io ^-zy 



x^C,2 + (1 ~ ^^) lii(l ~ a;)Li2(a;) + 5*1,2(3;) + x^lji^ix) 



-x{l — a;)Li2(a;) + j -x"^ -\ x ) In (1 — a;) 

1 



+-(l-a;2) 111^(1- 



X] 



(243) 



^ [^ 2;^ In (1 — zy) In (1 — xy) 

ay / dz 

^0 



I- zy 



r 



1 



-X — -X — - In (1 — a; 
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Jo ^- zy 



Jo 



I- zy 



--(1 - x) ln(l -x) + — {l- x^) 111=^(1 - x) 
— I — X X Liolx) H — a; C2 H — x G 

4 2 4 y ^^ ^ 4 ^^ 2 

+^{l-x')[S^,2ix)-U,ix) 
+^{l-x^)[+ln{l-x)U2{x)] (244) 

2^1,2(3;) + 2xU3{x) + 2(1 - x) ln(l - a;)Li2(a;) 

+-(l-a;)ln-'^(l-a;) (245) 

o 

-1 + 2(1 - x) ln(l - a;)Li2(a;) + 25i,2(x) + 2xLi3(a;) 



1 - - I Infl -x)--(l--] In^d - x) 



+ ii--ji„(i-.) 2 

+ |(l-i)ln'(l-a;) 



X 



(246) 



Jo 



1- zy 



X (Cs - C2) - (1 - 2x)Li3(x) - (2(1 - x) ln(l - x) 



+ (2 - x)S'i,2(a;) + (1 - x) ln(l - x)Li2(x) 

\ 
H — (1 — x) In (1 — x) — xLi2(x) 



(247) 



Jo 



I- zy 



2xC3 + 2(1 - x) ln(l - x)Li2(x) + (2x - l)Li3(x) 



^d f\~ yz^\n{l -y)\n{l-xy) 
Jo 



I- zy 



+ (3 - 2x)5i,2(a;) + -(1 - x) 111=^(1 - x) 
o 

<2 - - + x% + (1 - a:') ln(l - x)Li2( 
+ I x^ - - j Li3(x) + f - - x" j 5i,2(x) 



(248) 



X 



X — 1 



3x 



ln(l — x) — Li2(x) In (1 — x) 



1 



+ -(l-x^)ln^(l- 



X 



(249) 



'<;,/',, £i!!(i--) Ml--.) 

Jo 



I- zy 



:i-x) 



2C2 ln(l - x) - - ln2(l - x) - Li2(x) 



2Li3(l - x) + ln(l - x)Li2(l - x) 

+ a:C2 + 2C3 



+-ln^(l -X 
6 ^ 



(250) 



',,/',, £!l"(i-'-) Ml- -y) 

Jo 



1- zy 



^ + ^'Cs + ^ (1 - x') ln(l - x) [C2 + Li2(x)] 

3 1 

— -(1 — x) ln(l — x) — - (3 — 2x — x^) Li2(x) 

+ (l-x2)5l,2(x) + i(^X+^)c2 



40 



+- (2a;-3 + a;2)ln2(l - 



X) 



X] 



Jo 



1 - zy 



12 ^ ' ^ 

5x 7a;2 / 11 17 4 .A , , 
^ h ^ X + -x^ ln(l - x) 

69V 12 36 gy^ 

11 a; x^ x^\ ^ , ^ 2o, 
^ \ \ Li2(a:) + -x-^Cs 

18 3 6 gy ^ ^ 3 ^ 

I I ri^ rf" rf'~^ 

,ii + ^ + ^ + ^)ln-fl_a;) 
36 6 12 IF 



(251) 



+ i(l-x3)|ln(l-x)[Li2(x) + C2] 
+i(l-x3){ +i 111=^(1 -x) I 



(252) 



^ 1 [^ 1 yz^lnil — z)\n{l — xy) 
ay I dz 

Jo 



1 - zy 



2xC3 - ^ + xCs + 



4 V x) 



ln(l — x) 



"^7 [^ 1 yz^\\i{l — z)\\i{\ — xy) 

ay / az 

Jo 



1 -2;y 



(1 - X) 111^(1 - X) + -{I - X) 111^(1 - x) 

-(l-a;)Li2(a:) + 2(l-a;)5i,2(a:) 
+ (l-a;)ln(l-a;)[Li2(a;) + C2] 



(253) 



X 



(l-x2)ln=^(l-x) + {l-x^)Si^2^ 



X) 



12 

i(l-a;2)ln(l-x)[C2 + Li2(a;)] 

3 5 in, ., , , 2a; + 3;% 
—X In 1 — x H C2 

4 36 18a;/ ^ ^ ^ 4 ''^ 



+- (-3 + 2a; + a;2)Li2(a;) 



(254) 



'^r'^^z_Hl-yz)Hl-xy) 
Jo 



I- zy 



X (Cs + C2) - (1 - a;) [Li3(a;) - S'i,2(a:) + Li2(a;)] 



\ f\ y^\n{l-yz)\n{l-xy) 
Jo 



I- zy 



+ -{x - 1) 111^(1 -x) + -(1 - x) ln^(l - x) 
+ (1 — x) ln(l — a;)Li2(a;) 

a;C3-^ + ^(l--jln(l-a:) + -a:C2 



(255) 



+ (l-a;) 



ln(l - a;)Li2(a;) + - ln^(l - x) 
6 



3a; 1 



■1 + ^ + 



Ax 



In^fl -x) 



( \ 3a;\ ^ , , 

+ (l-a;) [S'i,2(a:) - Li3(a;)] 
41 



(256) 



Jo 



[I - zyf 



X (C2 - Cs) + a; ln(l - a;)Li2(a;) 
+x [5*1,2(2;) - 2Li3(a;)] + a;Li2(a;) 
— (1 — x) \v?{l — x) — xC,2 ln(l — x) 



(257) 



dy I dz 
Jo 



z^ hi{y) ln(l — xy) 
(1 - zyf 



3 / X \ 

-a; + a;^ln(l - x) [Li2(x) - C2] +a: f- + ij Li2(x) 

+x' [Si,2{x) - 2hh{x)] + ^ ln(l -x)- x\. 



X) 

V /"V 2;^ ln(w) ln(l - xi/) 5x 4x^ o^ srr^ / x ^.t • / m 

dy \ dz ,^^ _,, ^' = — + ^--x^C3 + x^[Si,2ix)-2U^{x)] 



(25J 







(l-zyy 



-G-r'-H'"'^-^>+(T-|-^')«= 



+ [ a;^ + -x^ + -X j Li2(a;) + a;^ ln(l - x)Li2(x) 

-C2a:^ ln(l - x) - - (l - x^) ln2(l 
6 



a; 



(259) 



io 



;i - zyy 



-2x (Cs - C2) - (1 - 2x) ln(l - x) [U2{x) - (2] 

-2(1 - a;)S'i,2(a;) + (1 - 4a;)Li3(a;) + 2a;Li2(x) 
-d-xlln^fl -x) 



(260) 



The following integrals were also required. They are understood in the sense of an analytic 
continuation and expanding in e to the order needed in the present calculation. 



/•I u 

/ dy y "^ ^ \n{l + uy) = \- u - {1 + u)\n{l + u) 

Jo ^ 



dyy '\n{l + uy) 



dyy* ^ln{l +uy) 
Jo 



+e [MLi2(-M) + (1 + m) ln(l + u) - u] 



2e 4 2 2 

,,2 Q „,2 

+e 



1-M^)ln(l + M) 



Y + i« - yLi2(-tx) + ^(1 - u^) ln(l + u) 



(261) 



(262) 



■H 2 Si 

\ \ (l + un n(l + u) 

- St • / X ""^ 4 2 5 

3 ^^ ^ 27 9 36 



/ dyy~^-'ln\l+uy) 
Jo 

/ rfi/|/-3-Mn2(l + M|/) 



i(l+M3)ln(l + M) 
9 



-2MLi2(-M) - (1 + m) ln^(l + u) 



+ -u^ - u{l + u) ln(l + u)+ u^U2{-u) 

c. ^ 



(263) 
(264) 



(265) 
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/■I u^ 1 1 

/ dyy'^-'ln^n + uy) = (l + u)(l-u + u^)ln'^(l + u)--(2 + 7u)u'^ 

Jo e 3 6 



u 2 

--(1 + u)(l - 3u) ln(l +u)- -u'^U2{-u) (266) 



t/y y Mn(l + Miy)ln(l + M2y) 



13 1 

■^MiM2 + TM1M2 - T ln(l + Ml) ln(l + U2) 



+- 



lit 



Li2(-U2) - Li2 
+ ln(l + M2)ln 



Ui 



+ LU 



U1 — U2, 

U2{l+Ul) 
U2 — Ml 



^1(1+^2) 
Ml — M2 



U2 



-y (1 + Ml) ln(l + Ml) + {Mi ^ M2} 



/■I 111 

/ c/y |/~^"^ln(l + Miy)ln(l + M2y) = — M1M2- -M2M1 - -ln(l + Mi)ln(l + M2) 
Jo 25 6 6 



--^ [m2 - (1 + M2) ln(l + M2)] - ^Li2(-M2) 



Ml 

' 6 

,3 



Mn 



M2 + (1-M^)ln(l+M2) + Y 



Ml 



H — ^Li2 

3 \Mi — M2 



_ m| / mi(1 + M2; 
3 V ""i ~ ""2 



3 V ^2 - Ml / 



(267) 



(268) 



X 



ldyy(l-y)U,^ (l-x)-,(l-y) 



— 7 — ^-TTT + - ln(l -x) ln^(l - x) 

9(l-x)2 9 ^ ^ 3 ^ ' 

1, , ,, ,, , 1 x(5x2 - 12x + 3) 

+- In X In 1 -X + -^— ^ 

3 9 1 — xr 



ln(x) . 
(269) 



B.5 Results for the Feynman integrals 



We will now present the results for all of the integrals appearing in equations (I106P to (I113p . 
We give the results up to 0(e°) in x-space. Logarithms, polylogarithms and Nielsen functions 
appear repeatedly in the expressions, so in order to save space, we use the following shorthand 
notation: 



Li = ln(x). 



L2 = ln(l-x), L3 = ln(l + x), Pi = ln(x) 



P2 = Li2(l-x), P3 = Li2(-x) 
P3 = Li3(-x), P4 = 81,2(3;), 

Let us define for X = A,B,C,E 



Pi = Li3(x), P2 = Li3(l-x), 

-R5 = 81,2(1 - x), Re = 8i,2(-x) . 



X, 



a,b 

1/IV2U3U4,U5 



dx x"(m2)4-^i2345+e ^^ . pY+bx^ 



a+b -^afi 

U11/2U3U4U5 



and 



:pa,b 
1^11^21^31/41^5 



dx (x" - 1) (m2)^-^^2^«+" (A ■ p) 



a+b fpafi 



Uiiy2U3U4,U5 



(270) 
(271) 
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The A-tjpe integrals are 



4' 



0,n 
01111 






4Pi-4C2--L? 



A, 



l,n 
01111 



2 1 

— [l-x+{l+x)Li] -- 



^C2Li - 6C3 - ^Ll - 2i?i + 8i?4 



a;) ( 4Pi - Li - ^L? - 4(2 



+ (l-x)(2-4L2) 



1 



+ ^(3a;-l)L2 + (l+a; 



4' 

■ri., 



0,n 
02111 



+2C2 - i^i - 2i?i + 8i?4 - — i:f 
P2 + ^i^? + C2 



+ (l-x) 



A' 



0,?i 
10111 



--(1 -x)- -(1 - x) [L2 -l] + {l-x) (3C2 + I6L2 - 8Ll - 16) 



^(1 - x) + - (1 - a; + xLi) + 8(1 + x) (P3 + L1L3 



A' 



0,n 
11011 



+ (3 + 5x)C2 + 4:xLl - IQxLi - 16(1 - x 
(-1) 



4^1 + ^Ll + 4Pi + 8P3 + 2C3 - 4L1P3 - 2L1P1 

5^ 2£ 



A 



l,n 
11011 



12^' + 2^^^^ 



-1)" I - 4 [^^1 + 1 - a;] - - 

I,. „ ..o . . „ 1 



2a; - 2 + -xLJ + (1 - 3a;)Li 
1 



, -(1 + 9x)Li + 4a;LiP3 - -(1 - a;)C2 + 2a;LiPi - -a;C2^i 

7 
xLl - 2^-iX - AxRi - 8a;P3 + (1 - 3x)Li + 2a; - 2 



A 
i' 
A 



0,n 
12011 

0,n 
11101 

l,rt 
11101 



(-1)" 

I0,n 
^11011 



^L? - Pi - 2P3 + 2-^ - L1L2 - 2L1L3 

2 1 + X 



— [1 - a; + Li 



£" 



2 - 2a; - -L^ + (3a; - l)Li 



+ 2&i. 



T0,n 



-i(l + 9x)Ll + i(l - a;)C2 + 2(3 + ^L^ - AL,P^ - 2L,P, 



+4Pi + 8P3 + (3x - l)Li + 2 - 2x 



O.n 



A 



12101 

0,n 
11110 

l,n 
11110 



^12011 



6C3 + ^C2i^i - ^Ll - 2Pi + 8P4 

-4(1 - a;) - -[(1 + a;)Li + 4(1 - a;)L2 - 2 + 2a;] + i(l - 3a;)L2 



4Pi-4C2--L? 



(272) 



ic2^i - 6C3 - 2Pi 
4L2-4L2-^C2 + 2 


(273) 






(274) 



(275) 



(276) 



(277) 

(278) 
(279) 



(280) 
(281) 
(282) 



:i - x) [All - 4L2 + 2) + (1 + a;) (Li + 2Pi) - ^H + 7a;)C2 



(283) 
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A 



0,n 
12110 



A. 



0,n 
22110 



P2 + -Ll + C2 

^4-^ + ^C2)^(l-.)-Q-l)D0(.) 



-22)i(x) 



1 — X 



TO,n 
^21110 



.l,n 
^21110 



/in,0 
^01111 



-- - - [2L2 - Ll] + 2P2 - 4L2 + 4L1L2 - -C2 



a: 



n,l 
01111 



-2R2 + 4i?5 

2 1 

— -\x + (l -x)L2] - - 



^- x\ C2- 4Ll- {1 + x)P2-^xLl 



4P2- 4(2-2^2 



1 
12' 



5" 



e 



-C2L2 — 2(^3 — -L2P2 + 77;'L2 



1 



i2-x)L2 + il-x)(^Ll + AC2-^P2 



+4:xLi - 2x 



-2x -il-x) 



+ ix - 4)Pi 



4L1L2 - i(2 - a;)L2 - 2xLl + ^(2 + (3x - 2)^2 
L2P1 + ^(2^2 - 2C3 + L,Ll + ^Ll + AR, - 2R2 



^10111 



7n,0 
^20111 



/<n,0 
^11011 



/.n,l 
^11011 



-- - - [L2 - 1] - 4 - 2Ll + AL2 - (2 



e^ e 



^ - - + 2 + \c^ 5{l -x)- -2)o(x) + 22)o(x) - 22)i(x) 



e^ e 



5 



7 



1 



-Ll + — Lt + 4i?i + 8P3 + 2C3 - 4L1P3 - 2L1P1 + —L\ + -C2L 

e^ 2e 12 2 

2 r 1 1 
— [xLi + 1 - x\ 



(3a; - l)Li - -xLJ + 2 - 2a; 



a;)C2 



+ IC2L1 - 2a;LiPi + 2(32; + 8a;P3 + (3a; - l)Li + 4a;Pi + —xLl 



— (1 + 9a;)Li - 4a;LiP3 + 2 - 2a; 



A 



n,0 
21011 



P2 - 2P3 - C2 + -L? + 2 — — 

■^ ^ 2 ^ 1 + a; 



2L1L3 



/>n,0 
^11101 



^n,l 
^11101 



. Jl _ 1 [412 - Ll - 2] - ^C2 - 7L? - ^i^Li - 2L2 - 2 



1 

A' 



1 — X 



X 1 + X 



-Ll - 2L2 + 1 



:i-x) 



lI + IlI + \c2 + i 



(284) 

(285) 
(286) 
(287) 

(288) 



(289) 
(290) 

(291) 
(292) 



(293) 
(294) 
(295) 



A 



n,0 
21101 



+ -(l + 3a;)Li 

4 + ^1 5(1 -x)- -2)o(x) - 22)i(x) + 1 + i^ I 1 



1 — X 



/in,0 
^11110 



e^ e 



(3L2 



1,2 1 

4 



L2 + 0C2 + Pi + 3L2 - 2 



Ll 

1 — X 



(296) 
(297) 
(298) 
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^11110 



Anfi 
^21110 



^ - - + 2 + ^C2 ) (^ll - x) - -Do(x) + 2Do(x) - 2Di(x) 

--Pi - --^2 ~ C2 



A 



n,0 
31110 



1 

2^2 



i - i + 1 - ^C2 ) <5(1 -x) - Do(x) - 2)i(x) 



-2+e 



a; 



n,l 
21110 



'3 1 

_2 ~ £ 

7 



2 + iC2). 



+ C2 + 2^2 + ^1 



£^ 



2e' y-A'^^ 



i) - i (Boll) - 1) + |l)o(i) 



A' 



0,n 
12111 



-2)i(x) - ^ - a;C2 - |i^^ + ^i^2 - xPr 

1 2 2 12 2 

b(\ -x)- -2)o(x) - -P2 - -^? - -C2 - -, ^ 

3e ^ ' 3 "^ ^ 3 3 ^ 3^ 3 1-x 



+ i(l-x)-^+^^L? + ^(l-x)-^+^^Li 



1 



i^i 



C2U(l-a;) 



\-\r 



4 2 , 

-P3 - - 1 
3 3 



2 



i^2 + 3C2 



-1 J 4 2 / 2 

^3 ' '"^3(l+x)2"^3(l+x) V(l + a;)2 1 



1 



X 



-iVi 



^? 



3(1 + a;) 



.n,0 
^21111 



i,^ + 3^2 + i ) 5(1 _ a;) + 2)i(a;) + f --^ - - ) C2 - ^-L\ 
2^2 8^ 2; ^ ^ ^ ^ V(l-a;)2 27^ 4 ^ 



;i-x)2 
1 



1 (L1L3 + P3) + 1 



1 



2(1 -x) 



P, 



U 



2(1 -x) 



{\-xf 



\x{x-2)L\^^^^U^\{x-2)UU 



(299) 



(300) 



(301) 



(302) 



(303) 



(304) 



The integrals ^inii; ^liiii ^^'^ ^liiii ^^^^ finite and appear multiplied by a factor of D — 4 = e 
in the expressions of diagrams [5^ and [5t. Therefore they only contribute to those diagrams 
starting at 0{e). For this reason we have chosen not to present them here, although they are 
needed in Eq. ( 1T24D . 

The P-type integrals are given by 



B 



0,n+l 
12011 



-1)" 



1 — X 



1 

2 



X. , , Q tXj X Ju 



-(5-x)Li + 2-2x 



+ 2x-2--(5-x)C2 



+P2--(7 + x)Lf + 2 



1 + x 



-Lx — ILxL'}, — 2P3 



B 



0,n+l 
12101 



B 



0,n+l 
12011 



(305) 
(306) 
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B. 



71+1,0 

21011 



5 , 7 + a; ^9 



^-^[(5-a:)Li + 5(l-x)] + ^C2- 
--(7 + llx)^Li + P2 - 2P3 - 2L1L3 - -(1 



X 



5, 



n+1,0 
21101 



-^-^(4xL,-xL, 






+ :t7^ r^i - -Ll 



B 



0,n+l 
12110 



[l-x) 
1-x 1 



X 

2' 



i<^ 



3a;) 



e' 



-{3x+l)Li + 2{l-x)L2 



i<^ 



X)C2 



1 



Rra+1,0 
-^21110 



1 - x)Ll - 3xLiL2 - 2xP2 - 2(1 - x)^^ + P^ 
1 3 7 1 .\ ... , X 1 



e2 2£ 



l + \C2]Sil-x) 



2^2 



e \ ^ ' 4 2 



3 3 7 15 3 3 a; 

+ -Do(a;) - Tiiix) a;C2 xL^ + -(1 + x)L2 xPi 

2 uv ; ^^ ^ 2 8 16 ^ 4^ ^42 



(307) 



(30^ 



(309) 



(310) 



The E-tjpe integrals are 



e: 



0,0 
01111 



4[^2-v ^ 



e 



3Ll - 3Ll - 2Pi + 2P2 



C2L1 + C2-^2 — ^-^1 



-5i?2 - 2L2L? - 2L1P2 - ^^1^2 + ^^2 - 3P5 - 7L2P2 + 3C3 

2 



£;, 



1,0 
01111 



\ [2xL2 - (1 + 2x)Li] - - 
-2xLiL2 + (4 + 2a;)C2 



1 



- 4(1 + x)Pi + - + 3a; Lt - 3xL 



+ (3a; + 4)LiLl - ^xL\L2 + (3a; - 2)C3 



-3a;P4 - xLxPx - [ a; + - J (2^1 + 2(a; + 4)^2 + (2 + 3a;)Pi 

-4(x + 2)L2Pi + —(1 - ^x)L\ + —xL\ - (8 + 3a;)C2i^2 
12 o 



i? 



0,1 
01111 



% [(1 - 2a;)L2 + 2xU\ - i 



+ ( ^ - 3a; ) L^ + 2a;C2 



4(1 - x)Pi + (4 - 2x)LiL2 + SxLJ 



x + 2) LIL2 - (x + 4)LiPi - xC2i^i 



9 



+ ( |x - 1 ) LiLl + {7x - l)L2Pi 



6x 1 C2^2 - -Ll - 2C3 



X 

2' 



l-22x 
12 



Ll + (3x + 4)Pi + (5x + 2)P2 



E, 



1,1 
01111 



2x 1 

— [(2x-l)L2-(2x + l)Ll]-- 



x 



2x(x + 2)C2 + ^(6x+1)L2 



X 



+2x(2 - x)LiL2 - 4x^Pi (6x - 1)P 



X , 



+ -(3 + 2x)C2i:i 



+x(2 + 3x)C3 - ^(5 - 2x)C2i.2 - ^(1 - 22x)Ll - x(l + 2x)LIL2 



(311) 



(312) 



(313) 



47 



^0,2 



-'^Ollll 



-|(8 + 5x)LiLl - x{l + x) [2Li + 71^] P2 + ^{l - 6x)Ll 
+5x{2 + x)R2 - 3a;(2 + a;)i?5 

^\x- 2x^Li + {2x^ -2x+ 1)^2] - - 
e^ e 



2 r2 



2a; + 4xLiL2 + Sx^Lf 



3x^ - X + - ) L^ + 2x(2 - x)Pi + {x- 2)L2 + {2x' - 4x + 4)P2 



-4a;Li - 4(2(1 - x) 



T i I \ 



- -C2 - ( 2a;^ + 3a; - - ) C2i^2 - 3a;^C2i^i + 2a; 



K 



2,0 
01111 



^^LlU + [ 1 - 2x + ^xM ULl + (3a;2 - 2x + l)PiL2 + x^PxL^ 

+4L1L2 - 4a;2L2P2 + 2a;2p2Li + 2xL\ + -^Ll - %^L\ 

+-^(22a;=^ - 2a; + \)L\ + (4 - 8a;)P5 + 3x2Pi + (Sx^ + 4x - 2)^2 

+ (4 - x)Pi + (2 - 3x)L2 + (4x - 2)C3 

2 1 

— \2x^L2 - (1 + 2x + 2x^)Li + X - 1] - - 



2x - 2 - 2x^LiL2 



-3x'^ - 4(1 + x)^Pi + - + X + 3x^ Lf + 4(1 - x)L2 + (1 + x)L 



+ (4 + 8x + 2x2)C2 
3 



2 

l + 7x 



- — C2 - X C2L2 - ( - + X + X ) (2^1 



1 — 3x 
+2x - 2 - -x^LlL2 + 2x^LiLl + 2x^L2Pi - x^LiPi + LJ 

-4(1 - x)Ll +[^ + ^-^^11 + ^x'Ll - (8 + 16x + 5x^)^4 
+ (3x2 ^ 4^ ^ 2)Pi + (1 + x)(Li + 2Pi) + 4(1 - x)L2 



+ (5x2 + 12x + 6)C3 
eZu = ^L2--{AC2-6Ll)+8Cs + ^Ll-7C2L2 



+{-ir 



- 4^1 - - (4P3 + 4L1L3 + L^ + 2C2) - III - 16P6 



^. 



1,0 
10111 



-I6L3P3 - SLiLl - 8C2L3 + 8C3 + 8P3 - 4L?L3 - 8L1P3 - C2i^i 

- (xL2 + 1 - x) - i [x (4C2 - QLl) + 8(1 - x) (1 - L2)] + 8(30; 

14 
+ (1 - x) {8Ll - 3C2 - I6L2 + 16) + —xLl - 7xC2L2 

' 4 2 

+ (-1)" — (xLi + 1 - x) - - (4xLi - 4xLiL3 - 4xP3 - 2xC2 
[e^ € 

-xLl - 4x + 4) - (3 + 5x)C2 + C2XL1 + 8xC2^3 + 16xLi 

-8(1 + x)LiL3 + 8XL1L3 + 4XL3L1 + I6XL3P3 + 8XL1P3 - 4xLi 

2 , 
+-xLl - 8(1 + x)P3 - 8XP3 + lexPe - 8(33; + 16 - 16x 



(314) 



(315) 



(316) 



(317) 



(318) 
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E. 



0,1 
10111 



- {xL2 + 1) - - (2a;C2 - 2^2 - SxLJ + 2) + (2 - 7a:C2i^2 + 8(32: + 4 



14 

y 



+^xL3 + 2Ll - 4L2 + (-l)"a; 



4^i + -(2C2 + i.? + 4LiL3 + 4P3) 



+C2 (i^i + 8L3) - 8C3 + -Ll + ALJLs + 8L1 (Pa + Lj) + I6L3P3 



-8P3 + IQRe 



(319) 



pl,0 

-^11011 



2V- -, J{2-x)Li + l-x]-2 + 2x + {2x-4:)Li-^^-^Ll 



— -il-x 

e 



2 

e 



+(1 - x) ( P2 - 2P3 - -C2 - 2L1L3 ) + (-1: 
1 /5 



- - (Li + 1 - x) 



e V2 



- -L^ + (4 - 2a;)Li - 2a; + 2 - -(2^^! - 2(3 - —L{ - 



7 



a: + 4 



12 



Ll 



^0,1 
-^11011 

-^11011 



+ {2x - 4)Li + 2L1P1 -{l-x)[ 2L1L3 + L2L1 + 2P3 + Pi + -C2 



+4L1P3 - 8P3 - 4Pi + 2a; - 2 

(-l)"^iion 
(l-(-inx 



1 /5 



1 



7 



(Ll + 1 - a;) - - ( -Lf + (4 - 2a;)Li - 2a; + 2 



2(2^^1 - 2C3 - —Ll + {2x - 4)Li + 2L1P1 + 4L1P3 - 8P3 - 4Pi 



X + 4 



p2,0 
-C'llOll 



X 



Lf - (1 - a;) -C2 + 2L1L3 + L1L2 + 2P3 + Pi + 2 



2 2 a; + 4 
- —(1 - a;) - - [(2 - a;)Li + 1 - a;] + 2(a; - 2)Li LJ 



+ (1 - a;) ( P2 - 2P3 - -C2 - 2L1L3 - 2 



+ (-1)" 



2 1 / 

— (2xLi + 1 - x^) + - 5xLl + (1 + 2a;) (1 - a;)Li 
e^ e \ 



-2(1 - x)M - 5xC2i^i + 2(1 - x^) + 2x^5 + lOxPi + I6XP3 

/I 1 

+ (1 - x) ( -(1 + 3x)C2 - a;P2 + 2xP3 + 2XL1L3 - ^(1 + '2x)Ll 

7 
+2C3X - (1 + 2x^ - 7x)Li + xLi (6P2 - 8P3) + 5xLlL2 + -xL\ 



6 



p0,2 
-'^11011 

pO,0 
-^11101 



2,0 



1 \n 771^, u 
-Lj -^11011 



1 — X 



1 I I 



-Ll --[-Li- 4L1L2 - 2P2 + -C2i^i + -:Li + 2LiL^ 



+2P2 - R5 



e \2 

(-1)" 
1 + x 



1 



12 
^Li - i (2C2 + 4P3 + SLJ + 4L1L3) 

-2C2L3 - C2L1 + 3C3 - 4L3P3 - 2L1P3 - hi - SLJLs + 2L1P1 



(320) 
(321) 



(322) 



(323) 
(324) 
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-2LiLl - 2i?3 - ARi - ARq 



E. 



1,0 
11101 



X 



1 — X 



%Li - - {]-Ll - AL1L2 - 2P2 ) + ^2^1 + ^Ll + 2LiLl 



e^ 



-2Ro — Rn 



e \2 

-11 



l + x 



12 



'^-^L, + -( 2C2X + l{x -^)Ll 



+AXL1L3 + AxPs 



I — X 



C2L1 + 2XC2L3 - (Sx + 2)C3 + 2L1P1 



^ ^-L\ + SxL^Lg + 2(3x + 2)LiP3 + 2x1^11 + AxL^P^ - 4Pi 



E' 



0,1 
11101 



12 
+AxRq - 2(3a; + 4)^3 

— 1^1 ^ ^ i 

e^ V 1 — X I e 



X 



l-x V 2 

X 



Lj - 2P2 - AL1L2 ]+2-AL2 + L 



C L 

+ — + — + 2LI + 2 + 

24 2l-xV2 



1 1 

-(2-^1 + 2L1L2 + T^-^i — R5 + 2-R2 



+i±^Li + (-ir ^ 



1 — x 



l + x 



^L^ + i (2C2 + 3Ll + 4P3 + 4L1L3) 



+C2i:i + 2(2^^3 + 2^1 + 3^i^3 + 2L1P3 + 2LiL^ - 2L1P1 + 4L3P3 



E 



1,1 

11101 



2P3 + 4Pi + 4P6 - 3C3 

2 



£^ 



+ 



2 / X ^ 

— { X + Li 



1 — X 



X 

e 



Li - AL2 



X 



2fl-a;) 



L\ - 8L1L2 - 4P2 



X 



1 — X 

l + 3a; 



-^2-^1 + TT^^l + 2L1L2 + 2P2 — R5 



X 



X 



+ -(2 + -jLi + 2xLi 



+x- 



1 — X 



Li + 2x 



.l)n 



l + X 



—x{l - x)Li - - 2x'C2 + 4^x'P 



5 — X 



X 



X 



+Ax^L^L3 - ^-^xL\ 1 + 2(1- ^)C2^i + (2 + 5a:)a:C3 + Y^(7 + ^)^\ 
-x^ (2(2^^3 + 3L2L3 + 2L^Ll + 4L3P3 + 4P6) - 2xLiPi + 4xPi 
+2a;(4 + 3a;)P3 - 2x(2 + 3a;)LiP3 



E- 



2,0 
11101 



X 



1 — X 

+2P2 ~ -R5 
+ (l-3a;)Li + 
+ (l-3a;)Li + 



|li - i \h.\ - AUL2 - 2P2 j + ic2i:i + ^i^? + 2LiL2 



-1)" 
,.2 



-^ 1-a; 



2 

l + a;2 
l + x 



12 



^1 -- 



S + x^ 



£ V 2 l+X 



-L\ 



^ (2C2 + 4P3 + 4L1L3) - 2(1 - x) j - i(l - X)C2 



i±^L?-2 
4 ^ 



4-3x2 7_a;2 2-5x2 
P3 — T-TTz ri^i -r—sz 



l + X 



12 l + X 



l + X 



(325) 



(326) 



(327) 



(328) 
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2-3a;2 



—— (4i?6 + 2(2^^3 + 2LiLl + AL^P;, + SL^Lg) + 2-- L^R 

1 + X ^ ' 1 + X 

^ ^ ^ :C2Li + ^^ (LiPi - 2i?i) - 2 + 2a; 



2(1 + a;) 



pl,0 
-^11110 



-'^11110 



1 + a; 
1 /. . 7 



9 



E- 



1,1 

11110 



- (L2 - ^i) - - ( 4C2 - -i.^ - 2Pi + 2L1L2 + -Li\ - -C2L2 + 7C3 

1 35 1 

--C2l>i - 2LiL2 - P1L2 + -L^ - P5 + 2Pi + -Lf - QR, 

2 1 / 7 \ 1 13 

^ (1 + i^2) - - ( 2 + 2C2 - 3L2 - -L2 - 2P2 ] - -C2 - y (2^:2 + SCa 

35 7 

+—Ll + 4L2P2 - 2P2 + i^i^2 + P1L2 -R5-P1- 3L2 + -LI + 2 

^"^ " '1 + 1) + 7 (2^2 + IlI + 3L2 + 4Pi - 6C2 -1^1-2 



— {L2 - L, 

—X 

-^-Ll - L2P1 + 3L2 + Pi + 2P2 - 2Pi + 8P4 + i?5 - 2 



13 
Y 



1 
2' 



1 
12' 



oC2^i + ^C2^2 + ^C2 - 9C3 - i^Ll - LiLl - AL2R2 - ^Ll 



35 

12' 



E- 



2,0 
11110 



2 1 
— {xL2 - xLi + 1 - x) 



2 - 2a; - 4a;Pi + 6a;C2 - 4(1 - a;)L2 



+ (1 - a;) 4Ll - AL2 - -C2 + 2 



-2xP2 - (1 + x)Li + |l2 - -xLl 

+ (1 + x) (2P2 - Li + 2L1L2) - 8a;P4 + 2a;Pi - a;P5 - 2a;P2 + X^x 

L? H xL^ H L? (2^1 + 3a;L2P2 a;C2^2 

4 ^12 ^ 12 ^ 2 2 



E, 



1,1 
21110 



-{V,{x) + 5{l-x)-l)-- 



3Ln 



X 



+l + (l + C2)5(l-a;) 



1 — a; 

7 



Li-3Di(a;)-2)o(a;) 



+ Y^^2 + ( - + 3a; ) C2 + (a; - 7)L 



-xRi - 2L2 - x^—^L,L2 + 2 + ("2 - ^C2 + 2C3 ) S{1 - x) 

7 



- ( 2 + -C2 ) 'Doix) + 2Di(x) + 7D2(a;) 



1,0 



-^11111 



1 — a; 



2a;P3 - 2P2 + 2a;LiL3 



X 



1 — a; 



Lj-i2- 3a;)C2 



+2(C2 + (3)^(1 -x) + 



(-1)" 

1 + a; 



- 2(a; + 2)P3 - 2- ^P2 



1 + a; 



-2(a; + 2)LiL3 - {x + 2)C2 + x-^LJ 

1 + a; 



nO,i 



-^11111 



1 — a; 



2(1 + 2a;) (L1L3 + P3) - (1 - 3a;)LiL2 - (1 - 4a;)Pi 



-L 



a;(5-7x) ^2 
■2(l-a;)2 



Lf + 2(C2 + C3)'5(l-:r) + 



1 + a; 



1 + X 



^ ^ 2a;(l - x) 
+2xL,Ls + ^ ^P2 + 2xP3 

1 + a; 



(329) 



(330) 



(331) 



(332) 



(333) 



(334) 



(335) 



(336) 
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E 



1,1 
11111 



1 — X 



x{l + x)LiL2 + x{l + 2x) (2P3 + 2L1L3 + Pi 



X I ox O I 2 

2ri-a;)2 



Lf - 4(22: + 2(C2 + (3)^(1 -x) + 



l + x 



-Ll 



2a;(l-a;) 
l + x 



x\A + x)^2 



l + x 



L{ + x{2 + x) (2P3 + 2L1L3 + C2 



E. 



2,0 
11111 



I — X 



+2x^P3 



x^ -2x- l)P2 + 2x^LiL3 + {Ax^ -2x- 1)C2 + 

-ir 



+ 2{C2 + (3)6(1 -x) + 



l + x 



+ {2x' + 2x - 1) (2P3 + C2 + 2L1L3) + (x' + 2x - 1 



a;(3 — 8x — 5x^] 
2ilTx) 



2 I O^ 1 \ p 



x{l + X 

2(1 -:r; 



-^? 



1 + a; 



(337) 



(338) 



The F-type integrals are given by 



0,0 



77lU,U 

-^01111 



71,0 
01111 

70,1 
01111 



(i-.)L-?^^-7G^^-*^^ + *«0-i«=^^ + 2«» 



1 

12' 



-— L^ + L2P2 + 2P2 - 4P5 



xP 



0,0 
01111 



4(^2 



-;l 



x 



I — X 



1 / 1 

- 4P2 - (4a;Li + (2 - a;)L2 - 2x) - 4(2 

£ V 1 — a; 



^(2^:2 - 2C3 + ^Ll + 4P5 - 2P2 - i^2P2 - ^P2 



+ -r^ ( ^^C2 + ^Ll + (2 - 3x)L2 + X [UL^ + 2L? + 2) 



71,1 

01111 
72,0 

01111 



70,2 
01111 



I — X 
•^-^01111 

2^0,0 
•^ -^01111 



2(l-x)Lo + 



x{2 — 3x) 
I — X 



;i-x)(iL^ + 4C2-4P2 



3x^ -lOx + Q 2a;(2 - 3a;) a;(3a; - 1) 

+ — — ^ L2 + ^ -Ll + ^ ' 



2fl-a;) 



I — X 



I — X 



3a;2-10a; + 6 . 9a;2 - 12a; + 5 ^ 9a;2 - 38x + 24 
-Li- 



[l-x) 
x{2 — 3a;) ^ ^2 1 



-Lo- 



2x 
I — X 

C2 



Ll 



I — X 
+ il-x) 



Ll + -L1L2 



2(1 -a;) " 4(1 -a;) 

3(a;2 - 6a; + 4) a;(l - 6a;) 



2fl-a;) 



-P. - 



2(1 -x) 



2C3 - ^C2^2 - ^i^2 + ^2^2 + 2P2 - 4P5 



pO,0 
-^10111 

-'^10111 

pO,l 
-f^ 10111 



1 — a; V e 
•^ -^10111 



4 - - (1 - ^2) + 2L2 _ 4L2 + C2 + 4 



I — X 



2-4a; 1 ,, , , 1 - 2a; , 
— — - - ((2 - 4x)L2 + 4a; - 1) ^(2 



(339) 
(340) 



(341) 

(342) 
(343) 



(344) 

(345) 
(346) 
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-^11011 



po,i 

-'^llOll 



;i - 2x)Ll + (1 - 4:x)L2 + 4:X- - 



X 



1 — X 



- ^Li - ^Ll - Ic^Li - 2C3 + 2LiPi - ^Ll 



2e 



+4.LiP^ - 8i?3 - 4i?i 



X 



1 — X 



Li 



12 



L^ + Li - 2 



2fl-a; 



1 — a; 



a; 



1 — a; 



2LiPi + 4L1P3 - 4Pi - 8P3 - ^(2^^! - 2C3 - ^Ll 



1, 1 + 9x^9 1 — 3x^ 
2' 4(1 - a;) ^ 1-x 



(347) 



(348) 



(349) 



-'^llOll 


~ ^-^11011 










(350) 


p2,0 
-^11011 




(351) 


p0,2 
-^11011 


8(1 - a;) 2(1 - a;) 1 - a; \ 


\ 
/ 




1 7 \ 1 — 3t 1 
+LIL2 - 6P1 - 8P3 - 2P5 - i^GLi --^Ll^+ ^^C2 - 2 - 3^ (352) 


pO,0 
-^11101 


1 


4--(i.i-4L2 + 2) + ic2 + 2L2 + lL? + i±^Li + 2 
e^e 2 41-x 


(353) 


l-x 


pl,0 
-'^lllOl 




(354) 


pO,l 
-^11101 


= -^(l-:^)^i°i?oi 


(355) 


-'^11101 


- rP°'^ 


(356) 


p2,0 
-^11101 


— -^ -^11101 


(357) 


pl,0 
-^11110 


X 


/ 2 1 17 \ 

- - - (2 - 3L2) - -C2 + -:Ll - 3L2 - Pi + 2 
\e'' e 2 4 / 


(358) 


l-x 


-^11110 


- 2^11110 2(1 - a;) ^' 


(359) 


-^11110 




(360) 


pO,l 
-^11110 


1 pi,o 1 . 
- 2^11110 2(1 -x)' 


(361) 


pl,0 
-^21110 


= -a;(l-x)-2+^ 


[2 /C2 ^ 1 

-- + 2- 1 + 2 5 
. ^ V2 / . 


T / 1 \ 

1 ( r 1 p 1 r^ \ 


(362) 


1-x V 2 y 


-f^lllll 


X 


- hi + 2^^ (L1L3 + P3) + -^ ( (3x - 1)C2 - xL,L 
2 L — X 1 — x\ 




l-x 


2 




~xLJ + {l-2x)P2) 
2 / . 




(363) 


-^11111 




(364) 


-f^lllll 


— ^-^11111 


(365) 


p2,0 
-^11111 












(366) 
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